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20.  Abstract  (Continued) 

surface.  Terms  beyond  first  order  (second  and  third)  are  included  in  order  to  explain  cross- 
polari/.ation  effects,  which  have  been  observed  experimentally.  Limited  results  on  fourth-order  terms 
are  also  obtained.  However,  when  the  analysis  is  carried  only  to  fourth  order  and  these  terms  are 
averaged,  they  do  not  contribute  to  the  cross-polarization  and  their  effect  is  seen  only  if  specular 
returns  are  present. 
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AN  ANALYTICAL  STUDY  OF  RADAR  RETURNS  IN  THE 
PRESENCE  OF  A  ROUGH  SEA  SURFACE 


INTRODUCTION 

The  analytical  study  described  in  this  report  was  originally  undertaken  in  an  attempt  at 
a  theoretical  explanation  of  certain  experimental  results  on  radar  return  from  ship  targets 
corrupted  by  sea  return.  In  the  process  of  carrying  out  the  study,  the  author  first  reviewed 
the  literature  with  a  view  toward  applying  the  existing  theoretical  results  directly  to  the 
problem  at  hand.  Previous  theoretical  analyses  by  S.  O.  Rice  at  the  Bell  Telephone  Labora¬ 
tories  (BTL)  [1] ,  a  group  of  investigators  at  the  Naval  Research  Laboratory  (NRL)  [2 
through  21] ,  other  investigators  in  the  Soviet  Union  [22  through  30] ,  and  still  others  in 
the  U.S.  and  elsewhere  [31  through  45]  had  produced  results  sufficient  to  provide  good 
estimates  of  equivalent  radar  cross  sections  of  the  sea  surface  that  explained  a  number  of 
experimental  observations  [3,46,47] . 

The  perturbation  approach  originally  introduced  in  the  radiowave  context  by  S.  O.. 
Rice  in  a  1951  paper  [1]  was  used  in  many  of  the  subsequent  analytical  studies.  It  was 
found  that  the  first-order  perturbation  results  did  not  explain  cross-polarization  effects  in 
the  return.*  These  effects  are  present  in  second-order  results,  which  were  obtained  formally 
by  Rice  [1]  and  later  by  Valenzuela  [6,20,21] . 

In  evaluating  the  problem  at  hand,  the  author  concluded  that  (if  the  perturbation 
approach  were  to  be  used)  second-order  fields  at  the  very  least  would  be  needed  to  explain 
the  experimental  observations  because  cross  polarization  effects  are  significant  in  this 
problem.  It  was  also  concluded  that,  in  the  presence  of  specular  returns  and  ship  target 
returns,  it  might  be  important  to  have  third-  and  fourth-order  fields  in  some  cases.  The 
reasons  for  these  conclusions  are  indicated  below. 

The  quantities  calculated  in  this  study  are  ensemble  averages  of:  (a)  twofold  products 
of  field  or  received  signal  voltages  (e.g.,  mean  radar  cross  section,  cross  correlation  function 
of  two  received  signal  voltages  corresponding  to  different  polarizations  or  different  delays, 
etc.)  and  (b)  fourfold  products  of  fields  or  received  signal  voltages  (e.g.,  cross  correlation 
of  two  received  signal  powers  corresponding  to  different  polarizations  or  delays,  etc.,  or 
the  power  fluctuation).  We  designate  terms  involving  products  of  fields  of  orders  j,  k,  fi, 
m  etc.  by  (jk9.m)  in  fourfold  products  and  (jjk)  in  twofold  products.  We  also  note  that 
under  the  assumption  of  Gaussian  statistics  of  the  surface  fluctuations,  odd-order  averages 
are  zero  and  specular  and  target  returns  are  considered  as  deterministic  and  of  zero-order 
in  the  perturbation  parameter.  Based  on  all  of  this,  the  twofold  product  averages  carried 


♦There  are  some  exceptions  to  this  statement;  e.g.,  Wright  [2,5]  using  first  order  theory,  considers  cross¬ 
polarization  due  to  surface  tilt.  That  effect,  which  is  also  covered  in  the  present  analysis,  is  a  type  of 
cross  polarization  from  a  different  mechanism  than  the  “quasi-crosspolarization”  discussed  later  in  this 
present  report. 
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to  fourth  order  will  contain  terms  (0000),  (0011),  (0002),  (0013),  (0112),  (0022),  (0004), 
and  (1111).  In  the  absence  of  zero-order  terms  (i.e.,  return  from  regions  with  no  specular 
or  target  return)  the  twofold  averages  are  confined  to  terms  (11),  (13)  and  (22)  and  the 
fourfold  averages  to  the  single  term  (1111).  In  this  case,  if  we  were  to  carry  out  the  twofold 
product  averages  to  second-order  and  the  fourfold  product  averages  to  fourth  order,  then 
the  first-order  fields  would  suffice.  However,  we  could  not  expect  cross  polarization  effects 
from  such  an  analysis.  These  effects  would  appear  in  the  twofold  product  averages  only  if 
we  were  to  carry  these  out  to  fourth  order  (i.e.,  if  we  include  the  (22)  and  (13)  terms)  and 
would  be  entirely  absent  from  the  fourfold  product  averages.  Their  appearance  in  the  latter 
type  of  averages  would  require  at  least  inclusion  of  terms  (1122)  and  (1113)  which  are 
sixth  order  and  obtainable  from  first,  second,  and  third  order  fields. 

If  zero-order  fields  are  present,  then  an  evaluation  of  twofold  product  averages  to 
second  order  would  include  (00),  (11),  and  (02).  The  (02)  term,  requiring  a  second-order 
field,  may  be  comparable  in  magnitude  to  the  (11)  term  in  the  part  of  the  return  that  is 
parallel-polarized  (i.e.,  the  VV  or  HH  return)  and  should  be  included.  The  HV  and  VH 
return  is  eliminated  in  the  averaging  process  for  the  (02)  term;  hence,  no  cross  polarization 
will  appear  if  we  truncate  the  series  beyond  second  order.  Again  the  VH  and  HV  returns 
average  to  zero  in  these  terms,  and  hence,  cross  polarization  does  not  appear.  In  the  case  of 
fourfold  product  averages,  all  three  of  the  terms  (0000),  (0002),  and  (0011)  must  be  included 
for  second-order  accuracy  in  the  VV  and  HH  returns. 

Since  the  second-order  fields  must  be  included  in  order  to  capture  the  cross¬ 
polarization  effects,  it  is  highly  desirable  to  use  them  to  calculate  the  fourth-order  terms 
(22)  and  (0022)  in  twofold  and  fourfold  product  averages  respectively.  If  this  were  done, 
however,  then  the  calculation  would  not  be  complete  without  the  inclusion  of  terms  (13) 
and  (04)  in  the  twofold  product  averages  and  the  terms  (0013)  and  (0004)  in  the  fourfold 
product  averages.  These  terms,  also  of  fourth  order  should  be  comparable  in  magnitude  to 
the  (22)  and  (0022)  terms  already  included.  Their  absence  would  give  an  incomplete  indica¬ 
tion  of  fourth-order  effects  in  the  VV  and  HH  returns  and  some  of  them  would  also  contain 
cross-polarization  components  (i.e.,  VH  and  HV). 

Now  let  us  proceed  to  a  discussion  of  the  analysis,  which  begins  with  an  electro¬ 
magnetic  wave  in  free-space  incident  on  a  sea  surface  with  two  scales  of  roughness.  There  is 
small-scale  roughness  (“ripple”  or  “capillary  waves”)  and  large-scale  roughness  (“swell”  or 
“gravity  waves”).  The  true  picture  of  a  sea  surface,  of  course,  is  a  continuum  of  scales  of 
surface  roughness.  Examination  of  previous  work  [6,7,39]  demonstrates  that  precise 
analytical  treatment  of  this  entire  continuum  is  extremely  difficult.  Hence  a  model  account¬ 
ing  for  the  two  limiting  scales  is  the  most  general  model  for  which  numerical  solutions  can 
be  obtained  without  great  difficulty. 

The  large-scale  roughness,  where  the  local  radius  of  curvature  of  the  surface  is  assumed 
to  be  large  in  terms  of  wavelength,  is  assumed  to  obey  the  “Kirchoff  approximation,”  i.e., 
in  the  absence  of  ripple  fluctuations,  return  signals  would  only  be  seen  from  regions  of 
normal  incidence,  just  as  if  the  mean  surface  were  an  infinite  flat  plane.*  Such  an  approach 


•As  already  indicated,  many  previous  theoretical  studies  have  made  use  of  this  perturbation  approach 
or  its  equivalent. 
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depends  for  its  validity  on  the  supposition  that  the  mean  surface  is  nearly  flat  over  many 
wavelengths. 

The  incident  and  scattered  waves  in  both  media  are  expressed  as  superposition  of  plane 
waves.  The  boundary  conditions  are  then  written  and  the  perturbation  theory  is  formulated. 
The  electric  and  magnetic  fields  are  expressed  as  superpositions  of  plane  wave  fields  at  an 
arbitrary  point  in  the  free  space  medium  and  are  then  determined  to  fourth  order  in  the 
perturbation  parameter.  The  electric  field  vector  at  the  radar  location  resulting  from  illumina¬ 
tion  of,  and  consequently  scattering  from,  a  small  region  on  the  surface  is  thereby  deter¬ 
mined  as  a  superposition  of  plane  wave  fields.  It  is  recognized  that  the  radar  is  well  into  the 
far  zone  of  the  illuminated  region  (denoted  by  S)  and,  hence,  the  backst  attered  radiation 
looks  like  a  plane  wave  at  the  radar.  Consequently,  a  single  plane  wave  is  extracted  from  the 
solution,  i.e.,  that  wave  propagating  from  the  region  S  (subtending  a  sufficiently  small  angle 
at  the  radar  to  be  equivalent  to  a  point  scatterer).  The  fie5-.!  vectors  from  all  such  regions 
illuminated  by  the  antenna  pattern  are  then  superposed,  i.e.,  the  solution  is  integrated  over 
a  set  of  spherical  angular  coordinates  ( 9,<t> ),  the  integrals  being  weighted  by  the  complex 
antenna  pattern  function  and  the  product  of  spherical  wave  Green’s  functions  e2^r/r2. 

This  approach  accounts  for  low  orders  of  local  multiple  scattering  due  to  ripple 
fluctuations.  This  is  inherent  in  the  perturbation  theory.  It  does  not  account  for  large- 
scale  multiple  scattering,  i.e.,  the  return  from  the  region  S  is  assumed  to  be  entirely  inde¬ 
pendent  of  the  illumination  of  other  regions.  The  latter  requirement  of  the  theory  would  be 
useful  for  treating  a  sea  surface  with  a  great  deal  of  large  scale  swell.  If  the  swell  is  small,  it 
should  not  be  a  very  significant  effect.  If  the  slope  of  the  swell  surface  were  very  small,  then 
illumination  of  any  region  by  secondary  radiation  from  other  regions  would  result  from 
wide  angle  scattering,  which  should  be  small  relative  to  direct  illumination  from  the  radar 
transmitter.  A  large  swell  surface,  however,  might  result  in  situations  where  a  given  region 
S  sees  specular  reflection  from  a  number  of  adjacent  waves.  In  this  case,  the  secondary 
illumination  could  conceivably  approach  the  same  order  of  magnitude  as  the  direct  illumina¬ 
tion.  Such  situations  should  be  rare  except  with  an  unusually  rough  sea  surface  and  hence 
the  neglect  of  this  effect  is  probably  not  a  serious  weakness  of  the  theory  in  most  cases. 

The  author  investigated  an  alternative  approach  used  by  some  workers  [39]  using  the 
“Stratton-Chu”  or  “Kirchoff-Huygens”  integral  equation  form  expressing  the  fields  at  a 
point  in  space  in  terms  of  the  fields  on  the  surface.  It  was  found  that  this  approach  gives 
answers  equivalent  to  those  obtained  by  the  direct  approach  outlined  above  once  the 
standard  approximations  are  invoked.  These  integral  equation  forms  are  merely  different 
ways  of  expressing  the  field  equations  and  not  solutions  of  these  equations.  It  is  still  neces¬ 
sary  to  solve  a  boundary  value  problem  in  order  to  determine  the  fields  on  the  surface. 

Hence  it  is  more  efficient  to  solve  the  boundary  value  problem  directly  to  find  the  field 
at  a  point  in  space  due  to  illumination  of  the  surface  by  the  incident  wave.  This  was  the 
approach  used  by  Rice  in  his  classic  1951  paper  [1] ,  and  the  use  of  the  integral  equation 
form  does  not  really  improve  on  that  approach  under  the  assumptions  made  in  our  model. 
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The  approach  used  here  is  essentially  that  of  Rice  [1]  which  has  been  used  by  sub¬ 
sequent  workers  (e.g.  Valenzuela  [4,6,20,21] ).  The  differences  between  the  present  develop¬ 
ment  and  that  of  Rice  are  as  follows: 

(1)  The  Fourier  integral  is  used  directly  instead  of  Rice’s  Fourier  series  with  a  later 
approach  to  the  continuous  limit.  In  terms  of  results  these  approaches  are  identical, 
although  there  are  questions  of  existence  of  Fourier  integrals  of  random  processes.  As  long 
as  one  understands  and  accounts  for  the  specialized  meaning  of  the  Fourier  transform  of  the 
random  process  while  carrying  out  the  development,  this  is  not  an  important  distinction 
between  the  development  in  this  report  and  that  originally  used  by  Rice. 

(2)  The  wave  incident  on  the  region  S  is  assumed  to  approximate  a  plane  wave  with 
arbitrary  polarization.  H  and  V  solutions  from  Rice’s  work  could  be  weighted  and  super¬ 
posed  to  treat  the  case  of  arbitrary  incident  wave  polarization.  This  could  have  been  done 
in  the  present  work.  The  principal  reason  for  not  doing  it  this  way  is  that  it  was  desired 
to  express  the  final  results  in  a  coordinate  system  with  origin  at  the  radar  and  z-axis  along 
the  antenna  beam.  Accounting  for  the  large-scale  swell  necessitated  the  use  of  a  complicated 
transformation  between  the  coordinate  system  used  in  the  basic  boundary  value  problem 
and  that  in  which  the  final  results  were  to  be  expressed.  In  view  of  this  coordinate  trans¬ 
formation,  there  was  not  significant  labor-saving  in  treating  horizontally  and  vertically 
polarized  incident  fields  separately. 

(3)  The  fields  are  calculated  to  fourth  order  in  the  perturbation  parameter.  Most 
previous  work  confined  itself  to  first-order  fields.  Rice  [1]  and  Valenzuela  [4,6]  calculated 
the  second-order  fields  and  used  them  to  evaluate  the  (2,2)  term  in  the  average  power 
(i.e.,  the  average  of  the  square  of  the  second-order  fields.  Valenzuela  calculated  second- 
order  doppler  spectra  in  later  papers  [20,21] .  The  reasons  for  our  inclusion  of  higher  order 
fields  were  delineated  earlier  in  this  section. 

Further  extensions  from  most  previous  theoretical  studies  are  as  follows: 

(1)  The  present  work  includes  the  effect  of  pulsing  or  other  forms  of  modulation  of 
the  transmitted  RF  wave.  Previous  treatments  are  generally  confined  to  pure  CW  trans¬ 
mission. 

(2)  The  present  work  includes  integration  over  all  incidence  angles  weighted  by  the 
antenna  pattern  shaping  function. 

(3)  The  present  work  includes  calculation  of  crosscovariance  functions  of  two  (pre¬ 
rectification)  received  signal  voltages  and  crosscovariance  functions  of  deviations  of  two 
received  signal  powers  from  their  mean  values  (where  the  two  returns  might  correspond, 
for  example,  to  two  different  polarizations  or  two  different  radar-target  geometries).  As 
degenerate  cases,  the  two  voltages  or  powers  may  be  the  same  but  with  different  delays, 
in  which  case  these  expressions  degenerate  into  autocovariance  functions  of  received  signal 
voltage  or  power.  With  zero  delay  difference,  the  latter  expressions  degenerate  into  mean 
power  and  “power  fluctuation”  respectively.  All  of  these  possibilities  are  contained  in  the 
averages  evaluated  in  this  work. 
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(4)  The  present  work  includes  effects  of  target  signals  in  the  ocean  environment,  i.e., 
it  considers  the  entire  “system”  problem  and  not  merely  the  CW  return  at  a  single  angle  as 
do  most  of  the  previous  theoretical  studies  of  sea  scatter.  The  target  return  and  sea  return 
can  be  considered  as  statistically  independent;  hence,  in  calculating  return  signal  power 
using  only  first-order  fields,  sea  return  and  target  return  add  incoherently  and  there  is  no 
coupling  between  the  two.  In  accounting  for  higher  order  terms  or  in  calculating  fourfold 
product  averages,  terms  arise  that  are  products  of  sea  return  and  target  return.  This  sort  of 
coupling  is  accounted  for  in  the  present  work. 

To  summarize  these  points  of  difference  between  this  work  and  previous  work,  the 
latter  has  stressed  the  calculation  of  the  scattering  cross  section  of  a  patch  of  slightly  rough 
sea  surface.  The  present  work  should  be  considered  as  a  “systems  study”  in  which  the 
traditional  theory  with  small  modifications  is  applied  to  a  system  consisting  of  (a)  a  radar 
whose  beam  illuminates  a  region  of  the  sea  surface,  possibly  a  large  region  and  (b)  a  target 
which  the  radar  may  be  trying  to  detect  or  locate.  The  analysis  accounts  for  possible 
pulsing  or  other  modulation  of  the  radar  signal,  the  shape  of  the  antenna  pattern,  and  the 
polarization  of  the  radar  signal.  A  certain  form  of  correlation  processing  is  done  on  the 
received  signal.  The  analysis  accounts  for  this  processing  and  focuses  on  calculation  of 
certain  statistical  averages  both  with  and  without  the  processing  and  both  with  and  with¬ 
out  the  presence  of  a  target. 

Previous  theoretical  analyses  of  scattering  from  a  rough  surface  have  been  examined. 
The  author  has  concluded  that  much  of  the  work  reported  in  the  literature  has  as  its  basis 
the  classical  paper  written  by  S.  O.  Rice  in  1951  [1] .  If  one  were  going  to  use  the  perturba¬ 
tion  theory  approach  to  the  problem,  it  would  seem  to  be  very  difficult  to  improve  in  a 
major  way  on  Rice’s  basic  work.  Subsequent  workers  have  usually  made  use  of  methods  of 
approach  equivalent  to  that  used  by  Rice.  The  present  author  has  also  followed  Rice’s 
approach,  with  a  few  small  extensions  and  variations  which  were  described  earlier.  Instead 
of  using  Rice’s  results  directly  or  those  of  Valenzuela  [4,6] ,  the  author  has  chosen  to 
rederive  the  results.  This  was  a  very  easy  (although  exceedingly  tedious)  task  and  was  for  the 
purpose  of  adapting  the  analysis  to  easy  extension  and  casting  it  in  a  form  that  could  be 
readily  prepared  for  digital  computation.  It  would  have  been  difficult  to  use  Rice’s  or 
Valenzuela’s  published  results  directly  for  the  calculations  to  be  made  in  this  study. 

Unfortunately,  it  is  true  that  the  third-  find  fourth-order  field  perturbations  are  very 
cumbersome  and  may  devour  computer  time  at  a  prohibitively  high  rate.  They  are  too 
cumbersome  for  hand  calculation  w’ithin  a  reasonable  time,  although  such  calculation  is 
simple  in  principle.  If  these  computations  are  prohibitively  expensive,  then  a  great  deal  of 
significant  information  can  still  be  obtained  from  the  first-  and  second-order  fields  alone. 

In  the  author’s  opinion,  it  is  still  worthwhile  to  have  worked  out  and  reported  these 
higher  order  terms  because  it  may  become  important  and  economically  feasible  to  compute 
them  at  some  later  time.  This  would  be  particularly  true  in  cases  where  cross-polarization 
effects  are  very  important.  In  such  cases,  the  second-order  terms  are  the  lowest  order  terms 
that  capture  these  effects,  and  third-order  terms  would  be  required  for  any  refinement 
on  the  simplest  results.  All  of  this  was  discussed  early  in  this  section. 

An  effect  not  accounted  for  in  the  analysis  presented  in  this  report  is  that  of  shadow¬ 
ing.  This  effect  assumes  increasing  importance  as  the  grazing  angle  decreases  and  the  surface 
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roughness  increases  [25] .  Since  numerical  studies  at  near  grazing  incidence  were  not 
contemplated  by  the  author,  no  analysis  was  done  to  account  for  shadowing.  However,  if 
such  cases  are  to  be  studied  based  on  the  theory  presented  here,  it  will  be  necessary  to 
supplement  the  model  with  further  analysis  that  includes  this  effect. 


THE  BOUNDARY  VALUE  PROBLEM;  CALCULATION  OF  THE  FIELDS 

Consider  a  rough  boundary  between  homogeneous  medium  #1  (free  space)  and  an 
arbitrary  homogeneous  medium  called  medium  #2  (e.g.,  the  sea).  The  mean  boundary 
surface  can  be  considered  as  an  infinite  flat  plane.  The  coordinate  system  to  be  used  is  a 
right-handed  rectangular  system  whose  origin  is  on  the  mean  surface,  whose  (x,y)  plane  is 
the  mean  surface  itself,  and  whose  z-axis  is  directed  into  free-space.  The  geometry  is 
illustrated  in  Fig.  1,  where  the  “mean"  incident,  reflected,  and  transmitted  waves  are 
shown  in  the  standard  way  for  a  case  where  the  incident  wave  is  a  plane  wave  and  the 
fluctuations  of  the  surface  are  designated  by  a  function  z  =  z(x,y,t)  which  is  a  random  func¬ 
tion  of  the  three  variables  x,y  and  the  time  t  and  has  zero  mean.  The  incident  wave, 
although  shown  as  a  plane  wave  in  Fig.  1,  will  be  arbitrary  and  possibly  pulsed  in  the 
analysis  to  follow,  but  will  be  specified  as  plane  at  a  certain  point  in  the  analysis.  The 
constitutive  parameters  of  medium  #1  (free-space)  are  denoted  in  the  usual  way  by  (e0  ,p0) 
and  those  of  medium  #2  are  (ec  ji0). 


Denoting  incident,  reflected,  and  transmitted  wave  fields  by  subscripts  i,  r,  and  t 
r<  spectively,  we  will  represent  all  field  vectors  by  column  vectors  and  use  matrix  notation 
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throughout,  having  found  that  this  leads  to  a  much  more  efficient  and  compact  analysis  than 
longhand  methods.  The  electric  and  magnetic  field  vectors  are  as  follows:* 


E(£,z  ,t)  = 


Eix^) 

r 

t 

Eiy{^,z,t) 


oo 

=  J' dooe~jcotEi  {g,z, co) 


Hiy{z,z,t) 


Hiz(p,z,t) 


oo 


where  q  =  (xy  )  and 


K  <0 


K  <0 

yz  (}) 


*)  I,(4.w) 


where 


4  = 


Pz  -  U,  (0  =  +  y/  1  -  ft,2  -  fly2 

72  =  7,(0  =  +  v/f'2  -  0*2  -  Py2 


U<4,z,0))  =  f-jdQe1" 


(0 

^  (0 
72  (0 


r  r 


t  t 


(2.1a) 


(2.1.b) 


(2.1.c) 


*1  (£)!,(4,w)  (2.1.d) 


•Note  that  some  of  our  notation  is  patterned  after  that  of  Rice’s  classical  paper  [1], 
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(where  fl  =  shorthand  notation  for  (0  ,o>)) 


ffd%  j(P,.  (Q)g  (g)e'Jc 


+  Pr(Q)Er(n)e  ,  ] 


(2.4.b) 


where 


[Pf  (Q)g,  (Q)e  ~'c  V5)z(4'f)j|e;£  wtl  =o 


0  -1 
x  0  z, 

*y  ’**  0 


At  this  point  (following  Rice’s  treatment)  we  invoke  the  divergence  equation  for  the 
electric  field,  as  follows: 


V  •  £  =  0 

Using  the  electric  field  representations  (2.1.a,c),  Eq.  (2.5)  takes  the  form 


ffdV  [jf  j  git  (A)  +  Pygiy  (Q)] 


-m 

+pz  <g> 


s. 


[ -y*W 

Multiplying  (2.6)  by 


(§)}*  JUtexp  ($<4'%) 


-7Z(0) 


=0 


-4-  <£'  -e.  - «*> 


integrating  over  all  /)  and  t  and  invoking  the  relation 

oo 


If 


djo< 


-  (2tt)2  (^)2  S(0X  -  Px'my  -  0y') 


(2.5) 


(2.6) 


(2.7 .a) 


-OO 
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OO 

J" J'dt  w  )*  =  2j r  6(co  -  c*/) 


(2.7 .b) 


we  obtain 


£.  (12)  =  +  —|i3  E-  +  Q  E-  |  , 

atz  ,  n  rx  IX  My  iy> 

vz  r  r 


r  . 
t 


r2 

1 


(2.8) 


We  can  now  construct  a  new  matrix  relationship,  i.e.. 


fi  = 

r 


£ 

IX 

r 

t 

% 

r 


^2 

-  _I  |[fi  £.  +8  £•  1 

p  ]lHx  IX  Hy  »y* 


r  r 
f  / 


where 


(2.9) 


Af;  = 


1  0 
0  1 


1 

0 


0 

1 


UjL  M  =  ,u  1  Af=ux 

/M  /M  ’  r  /^\  ’  '  /f*\  MA 

w  Vz)  \  h)  \  h)  W  \yj 


1  0 
0  1 


Stm 

r 

t 


E- 

ix 

r 

J 

E- 

<y 

r 

t 


I 


‘■vs-wy* iara**»**!**»*“v ' 
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From  (2.1.a,b),  (2.4.a,b)  and  (2.9)  we  construct  the  following  matrices. 


NMi  = 


-  B  z 
Hx*y 


(-t) 


(>*¥) 


$xzx  j  ^yZx 

X 


-  Z„ 


(2.10.a) 


NMr  = 


NMt  = 


B  z 
yx  y 


(-$) 


6  z 

H*y 


(■•t) 


(■-¥)■ 


B  z 
Hy  x 

X 


-  z_ 


B  z 

x 


-  z 


(2.10.b) 


(2.10.C) 


At  this  point  we  adopt  a  notational  convention  which  will  be  used  many  times  throughout 
the  development  to  follow,  wherein  a  prime  will  be  used  to  indicate  an  interchange  of 
x  and  y  components,  e.g.,  if  a  =  0xzy ,  then  a1  =  ($yzx .  Continuing  the  construction  of 
needed  matrices  with  the  aid  of  this  notational  device,  we  have 


NPiMi  =  /  = 


'll 

'l2 

hi 

hi 

hi 

hi 

(2.10.d) 
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where 


We  can  now  use  (2.10.a, ...,  f)  to  construct  a  4  X  4  matrix  equation  consisting  of  the 
x  and  y  components  of  Eqs.  (2.4.a,b)  (with  the  factor 

1 

>/i+K|2 


removed)  as  follows: 


//« 


dnAiCljMEWte^  t(^*fi)~  CfJ 


#9  B(5,P,0  £x(«)c 


>r  l(g*£)-  c<\ 


where  S2  symbolizes  the  set  of  variables  (Px  J3y  ,cj)  and  where 


■A(Slj),t)  - 


11 

*12 

*13 

*14 

12 

*11 

-*14# 

~*13 

31 

*32 

*33 

*34 

i 

32 

*31 

*34 

*33 

"511 

512~ 

-fi12 

)  ~ 

»«4 

CO 

a> 

B32 

lS32 

531 

5l 


•ft] 


(2.11.a)' 


(2.11.b)' 


(2.11.c)' 
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(2.11.d)' 


(2.11.e)' 
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B 


11 


Vxzy  -%&zz 

h  € 


We  now  invoke  our  perturbation  theory,  in  which  terms  of  A  and  B  not  involving 
z,zx  or  zy  are  considered  zero  order.  Thus,  using  f  as  an  ordering  parameter,  we  have 

A(«,p,f)  =  A<°>(fl,w)  +  (2.12.a) 

— ►  — ►  — *•  — ►  — ► 

B(J2jP,0  =  B<°>( B)  +  fB*1  >(/W,t)  (2.12.b) 


A<°>(g,w)  = 


0 

-1 

0 

1 

1 

0 

-1 

0 

A'°> 

A31 

A(°) 

A32 

A(0) 

A33 

A(0> 

34 

A(0) 

32 

A(°r 

*31 

A(0> 

*34 

A(0>' 

33 

(2.12.a)" 


where 
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-B  B 

A(o)  m 

34 


B<°Hi 3)  = 


0 

-1 


1 

0 


A°)  AO) 

31  A32 

AO)  AO)’ 

32  31 


A{1) 

All 

A{1) 

A12 

*(1) 

A13 

A(1) 

^14 

-A{1)' 

A12 

AY 

AY 

-A(1)' 

^13 

A{1) 

A31 

48 

*8 

A(l) 

A34 

aM 

AiV 

AiV 

AiV 

f 


:r, 


4  . 


where 
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tj  /  y2  - 

4V  =-Vy^rV-(— 


2  _  /?2 


4V  =  0* 


2ye 


-j?y,2  0; 


-^4 


®xV  B,iV 

-■ft"  -■IV 
■ft’  « 
■ft"  ■ft" 


rfl)  =  ^z  c-'rM 
fl»  0, *» 


,(1)  = 


^y  -tfrP* 


12 


£<1) 

^31 


d(!) 

S32 


2  6 


A  y 


[-J-P  Z  Py  (  -i-P  z  \1 

v  ‘ 1 1(' e  *-1)J 


Expanding  ^  in  a  perturbation  series, 

|=|°)+f|(1)+r2|(2)+.... 


(2.12.H)' 


(2.13) 
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and  employing  (2.12.a,b)  and  (2.13)  in  2.11),  we  can  write  (2.11)  in  the  form 
ffdifd^e7  [£*P)  -  ct]  ^((M^IO)^  _ 

+  f[>l(0>(^,co)  f^fi.Gj)  -  (B<^fl,co,p,f)f(.(/),co)  -  ^1>(0,cu,p,O$O)(|3^))] 

*r  — ►  —*•  — ►  — ►  — ►  ■*>  — ► 

+  r2[A(0)(^co)f  2)(^>to)-  (-  A^^t)^1^))]  + ...  (214) 

+  rM(0)^)|(n)(|,CO)  -  (-  ^1)(|,w!j0,f)|(n  •  X)%03))]  +  0(?"  +  !>)  j 
Multiplication  of  (2.14)  by 


-7jr  I((J'.P)-  cf) 
1  c  -♦  -* 


and  integration  of  product  on  jO  and  t  from  to  +°°,  with  the  aid  of  (2.7.a,b)  yields  the 
following  set  of  perturbation  equations  (after  interchanging  (0',u/)  and  (P,cj)  and  changing 
order  of  integration)  * 


^(°)(|,w)f  °)(|,w)  =  B<°>(£)  |(£<u) 


oo  oo  oo 


a(°)(p,oj)e(i>(p,oJ)=  — - 


(2jt) 


-  oo  — oo  -oo  — oo 


(2.15.0) 

~  ojJ3)'jP~  c(u>'~  <jj  )f  ] 

(2.15.1) 


(^(gV^of^V)  -  ^(1)(£',cu',p,t>^0)(r,w')) 


oo  oo  oo  oo 


A(0)(£,u>)|(2)(£,gj)  =  — ^  J'dtecl  U^  ^ 

^  ^  OO  OO  oo  oo 


-c(w'- w)<] 

(2.15.2) 
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(-  i4(1>(gV;Jg,t)f(n‘  1)(^',co'));n=3,4, ... 
Solution  of  (2.15.0)  yields 


where 


R(0)(t 3,^)  =  IA<°>(/3,<o)]  -1  5<o><0,w), 
and  the  inverse  of  A(°)  is  given  by 

T  *11 


11 

^12 

b13 

*14 

> 

12 

-*11 

*14 

*11 

11 

*32 

b13 

*14 

i 

32 

-*11 

*14 

*13 

where 


bll  ~  Px^y (Pz~yz)> 

b12  =  VzWy2> 

b13  ~  ~  +  @y  )» 


614  =  My>and 

b32  "  *  T,  -  tf*  -  7*)0y 
Solution  of  (2.15.1)  follows  from  (2.16.0),  i.e., 


I1  ><fM  -jfjfa  *»,  <*  *,  . 

£( 1  >(0,00/^  .COj  ,f  j  )|(-  (JJj ), 
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where 


*(1  >(& .«&  .Wj  )  =  —  M<°>(£  CO)]  -1  j  B( 1  ><&  ,Wj  ^  ,tj ) 

(2tt)3 

From  (2.15.2),  (2.16.0)  and  2.16.1), 

J2)(g.«)  d“i  d£i  dti  dJ2  dtJ2  d(>2  dt2  - 

""!>•£! 'C(W1_")*1 wifii)*fi»"c("*-"i>f2l_  (2.16.2) 

* ( 2  )(|,<o/|  1  ,U!  /|2  ,W2  ;pj  .tj  /p2  ,t2)  (P2,w2 ), 


where 


*( 2  )(|,w/|1  ,w1  /| jj  .Wjgj  ,tx  /g2  2 )  = 

”  (2^®  [‘A<°)^,W)l~1‘A(1)(j|1’&,1;-51*t1)^(0)^1.Wl)l  1|B<1)(J2,w2iS2d2) 
-A(1)(|2,W2ig^2([A(0)(|2,w2))-1B<0)(|2^2>i- 

The  general  (nth)  term  of  £  is  (from  (2.1 5.n)) 

fn)(g>w)  «jffa  dwl  d4l  d'l  dJ2  dcj2  dM-2  dt2  -  di«  dw„  d4n  dtn  ~ 
Jtwi£i  “  “  e(«x  -  fcijtj  )'i2"  c(w2- 


wn-2^,-2)^n-l-CK-l-  Wn-2>Vl 
+  K&  “  "n-lln-l)^  ~  CK  “  W«-1>U  - 

n  ,OJi  /j32  ,a>2  /  —  /^„  >OJn’$1^1^2^2/  — 


(2.16.n) 
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where 


_  (_!)«  - —  [A<°)(flIcJ)r1A<1)(01,o;lip1,f1) 

-  3  n  -*  -*1  1  -*1  1 


(27T) 


[^°)(|1,oJl)r1^1)(fl2,W2;p2,t2)...  M<°>Cg2,co2)l  -1^<l)<g3,co35e3,r3> ... 


IA(°%-  2  2 )}' ^(  1  x 


The  vector  E  is  not  of  direct  use  to  us.  It  would  be  preferable  to  have  the  reflected 
wave  vector  Er.  To  obtain  Er  we  invoke  (2.9)  and  note  that 


E  =  ME  = 

-»r  r  ->r 


1 

0 


0 

1 


0X  0y 

V"V 


"l 

0 


0  0  0 
10  0 


O' 


'fx 


<y 


(2.17) 


where 


Q(g)  = 


0 

1 


1 
0 

0*  -V 

v-v 0 


0  0 
0  0 


The  unit  vectors  corresponding  to  horizontal  and  vertical  polarization  are  defined  in 
Appendix  I  (Eqs.  (I-24.a,b)  and  (1-27))  and  are  denoted  by  By  and  respectively.  The 
specialized  meaning  of  horizontal  and  vertical  polarization  in  this  report,  different  from 
the  standard  meaning,  are  explained  in  Appendix  I. 
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The  vertically  and  horizontally  polarized  components  of  §r  are  denoted  by 

ErV  =  £.V  Jr’  and  (2.18-a) 

ErH=fHMr-  (2.18.b) 

In  all  that  follows  we  will  adopt  the  notational  convention  introduced  in  Appendix  IV 
wherein  all  explicit  indications  of  dependence  on  co  are  dropped. 

We  now  relate  E i  to  EiV  and  EiH,  the  vertically  and  horizontally  polarized  components 
of  the  incident  wave  vector  From  (2.8)  and  (2.9) 

!,  (0  =  (0-  (2-19) 


where 


(|) 

<0 


L($) 


^Hy  ~^Vy 
fiVx 


and, 


o'  =  o 
kHx  kHx 

y  y 

Vx  Vx 

y  y 


*Vz 

+  'J~X' 

r z  v 


X 

y 


The  remaining  analysis,  detailed  in  Appendix  IV,  results  in  expressions  for  vertically 
and  horizontally  polarized  electric  field  components  at  the  receiving  point,  given  both 
vertically  and  horizontally  polarized  transmitted  fields. 


THE  RECEIVED  RADAR  SIGNAL 

We  denote  the  received  signal  at  the  radar  by  the  symbol  u.(t).  This  represents  a 
quantity  proportional  to  the  electric  field  component  along  the  direction  of  polarization 
selected  by  the  receiver.  The  subscripts  A  and  a  denote  received  and  transmitted  polariza¬ 
tions  respectively.  We  can  use  the  standard  representation  for  narrow  band  signals,  i.e., 

uA.MmRe\UAu<t*iU°t  f>  (3-1) 
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where  UAg(t)  is  the  complex  envelope  of  uAa(t),  cjq  is  the  central  frequency  of  the  signal 
spectrum  and  UAa(t)  varies  slowly  compared  with 


The  most  general  quantity  that  will  ultimately  be  needed  involves  the  product  of  two 
signals  possibly  with  relative  time  displacement  At  and  different  polarizations,  time- 
averaged  over  an  RF  cycle  (where  an  overline  indicates  the  time  averaging  process,  and  of 
course 


e 


*2  j0Jot 


0. 


This  is  given  by 

VAoBbi‘M)  - \ReiUAaBb(tAt)e’“-\ 


(3.2) 


where 


Va.bi.UM-Va'I'Ki.V**). 

and  where  B  and  b  represent  received  and  transmitted  polarizations  that  may  be  different 
from  A  and  a. 

A  special  case  of  (3.2)  is  the  signal  power,  or  equivalently  (if  proportionality  constants 
are  set  equal  to  unity)  the  quadra tically  rectified  signal  voltage, 

-\  *«Va.a.W  -  \  6a.a. ('.»>.  <3-3> 


where 


‘WwHtfA.wl*- 

Another  quantity  of  interest,  in  its  most  general  form,  is  the  product  of  two 
quadratically  rectified  signals  with  relative  time  displacement  At  and  possibly  different 
polarizations.  This  is  obtained  from  (3.3)  and  is  given  by 

VAaBbW)  =  UAaAa(t,0)UBbBb(t  *  At, 0)  =  |i;Aa(t)  I2  |t;Bfc(f  +  At)  I2,  (3.4) 

neglecting  tire  factor  1/4  which  does  not  affect  the  results  significantly. 

We  now  invoke  the  simplifying  assumption  that  the  incident  wave  is  a  GW  signal, 
amplitude  modulated  by  a  pulsing  waveform  of  arbitrary  shape,  but  whose  bandwidth  is 
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small  compared  with  RF.  More  accurately,  we  should  solve  the  boundary  value  problem  for 
each  frequency  co  contained  in  the  incident  waveform  and  integrate  on  to  over  the  entire 
spectrum.  This  process  can  be  approximated  for  practical  purposes  by  treating  the  signal 
as  if  the  incident  wave  were  CW  with  frequency  a>  =  oj0  and  multiplying  the  resulting 
scattered  CW  field  by  the  pulsing  function.  This  implies  that,  within  the  small  band  covered 
by  the  pulse,  the  scattered  field  waveform  is  approximately  independent  of  frequency. 

This  simplified  model  is  implemented  through  the  introduction  of  a  function  g(ra), 
which  weights  different  contributions  to  the  return  in  proportion  to  the  delay,  which  in 
turn  is  proportional  to  r0 ,  the  distance  from  the  radar  to  the  illuminated  point  on  the 
scattering  surface. 

Another  weighting  function  on  the  signal  is  the  antenna  pattern  function  A 03,03) 
which  weights  contributions  from  different  angular  positions  (03,0*3)*.  Still  another  is 
the  factor 


(c2'7ro/r2) 

present  in  any  two-way  radar  signal  return  from  a  point  target. 

Allowing  fo.  the  possibility  that  the  pulsing  function  and  antenna  pattern  function 
may  differ  for  different  polarizations  we  will  place  subscripts  Aa  on  them  and  construct  an 
overall  weighting  function 


FAa($  =  * 


['■O(03>*3’f>l2 


(3.5) 


where  the  radar-surface  geometry  fixes  the  functional  dependence  of  r0  on  the  angles 
(03,03)  where  ^<,(03 ,03)  is  the  one-way  field  pattern  of  the  antenna  for  the  polarizations 
denoted  by  Aa,  where  k  is  a  constant  containing  radar  parameters,  and  where  the 
possible  time  dependence  in  r0  arises  from  the  expression  (Appendix  III,  Eq.  III.  24), 


ro-?o  +  Aro’ 


(3.6) 


where  r0  is  the  distance  from  the  radar  to  a  point  on  the  horizontal  (calm)  sea  surface, 
given  (through  manipulation  of  (III.25)  in  Appendix  III)  by 

hR 

r  = - - - 

0  cos  0g  cos  7  -  sin  0g  sin  7  cos  0g  ’ 


*See  Appendix  I  (beginning  below  Eq.  1.14)  for  definition  of  the  angles  (dj.^g). 
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where  is  the  radar  altitude  and  7  the  beam  angle  relative  to  vertical,  and  where  (from 
(III.24)  or  (III.26)) 


A ro=-H(x'y)- 


(3.6.b)' 


H(x',y')  represents  the  mean  sea  surface  as  a  function  of  the  horizontal  coordinate  (x'y'), 
hR  is  radar  altitude  and  fQ  is  given  by  (3.6.a)'. 

If  H(x'y)  =£  0  in  (3.6.b)\  i.e.,  if  the  mean  sea  surface  is  not  perfectly  horizontal 
everywhere,  then,  since  (through  manipulation  of  (IH.lS.f)  in  Appendix  III) 


(sin  0g  cos  0 3  cos  7  +  cos  0 g  sin  7) 

X  ~  X  n  ^  h  p  "  f  •  nl  • 

"  (cos  0g  cos  7  -  sm  flg  sm  7  cos  ) 


(3.6.C)' 


hR  sin  0^  sin  ^ 

y  =  y  — - - f 

R  (cos  0g  cos  7- sin  flg  sin  7  cos  ^3) 


(3.6.d)' 


where  xR  and  y'R  are  the  horizontal  coordinates  of  the  radar  relative  to  a  fixed  origin  of 
coordinates. 

The  time  dependence  in  Ar0  arises  from  a  possible  time-dependence  in  x'R  and  y'R 
which  would  appear  in  x',  and  y'  through  (3.6.c,d)'  and  hence  would  appear  in  H(x',  y') 
and  in  (3.6.b)'.  This  would  occur  through  horizontal  motion  of  the  radar.  Vertical  radar 
motion  would  manifest  itself  in  time  dependence  of  h R ,  which  would  render  rQ  time- 
dependent  through  (3.6.d)'  and  would  also  appear  in  Arc  through  (3.6.b,c,d)'.  Another 
possible  source  of  time  dependence  is  a  time  variation  of  the  beam  angle  7,  which  would 
appear  in  rQ  through  Eqs.  (3.6)  and  (3.6.a,b,c,d)'.  If  the  mean-sea  surface  were  perfectly 
horizontal  everywhere  and  the  beam  direction  fixed,  then  r0  and  hence  FAa  would  be 
independent  of  time. 

We  now  specify  the  form  of  the  radar  return  signal  as  a  superposition  of  returns 
from  patches  at  position  (rQ,  6 3 ,  ^3 )  relative  to  the  radar,  each  return  weighted  by  the 
function  FAa(0 3,  <P\,t)  as  given  in  (3-6)  and  a  factor  SAa  derived  from  the  perturbation 
solution  of  the  boundary  value  problem.  The  results  of  this  line  of  reasoning  are 


JAao(t) 


n  ZTT 

=fdQ'z  sin0'y'd0'FAa(0',0',t)SAao(ao(0'4,f)),  (3.7.0) 
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it  itt 

UA* 1«>  -/ “f>3  *>  «3 / »'  FaJ9> '3^SX.|(S.  <93.^.<»  - 


(3.7.1 


2gj 

&  <*>**». 


it  tilt 

UA,2<‘>  ~f  «k  *  *i/ #k  FA^y>ffdKSA.2tSc  4) 


(3.7.2 


Z(k)Z(-— 4,(03,03,0-*), 


ff  2ff 

*/ *>  / *»k  ffdi2SAo3  •■ 


<&,  <«i.<,i,‘).'ii.i2)Z(fe1)Z(fe2)^Z( - 4,(«k^kA- ij  -  *2),and 


(3.7.3) 


VUw  > 

I  hi  -  c-  H 

i  I  '  ,  ^  ‘  *5;*l 
f  r  'fe'iSl 


«  zjt 

^1.4(0  =/ i03  sin  03 / ^Aa^'z^ffel  ffdk2  ffdHz  S, 


AoA  - 


(3.7.4) 


^So^3’^3’^’  "• 

z<-^  s.(«i>i^-S»-4.-4»>. 

where  a0  is  defined  in  Appendix  I,  Eq.  (I.27.a),  as  the  projection  on  the  x-y  plane  of  the 
unit  vector  directed  from  the  radar  to  the  illuminated  surface  point. 

Note  that  the  a0  in  (3.7.0, 1,2, 3,4)  is  indicated  as  a  function  of  the  angles  (6'3,<p'3)  and 
a  possiblp  function  of  time.  The  time-dependence  arises  again  through  possible  motion  of 
the  radar  and  beam  and  departures  of  the  mean  sea  surface  from  perfectly  horizontal.  To 
show  this,  we  invoke  (I.27.a)  and  (1.28),  dropping  subscripts  in  (1.28)  and  also  Eqs.  (I.8.a,b) 
and  (I.10.a,b)  with  the  results 

4>  =  (ao*’aoy)’  (3-8) 
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where 


<v  = 
ox 


— 7 - . — p  jsin  do  cos  <t>'o  [cos  4>  cos  7  -  I V'H  I  sin  7] 

v/l+l  V'/jr  (3.8.a)' 

-  sin  0g  sin  0g  sin  +  cos  0g  [cos  $  sin  7  +  |V'//  j  cos  7]  j ,  and 


aoy  =  j-  sin  0g  cos  <t> 3  sin  <J>  cos  7  sin  0g  sin  <j>'3  cos  4> 
-  cos  6 3  sin  4>  sin  7^. 


(3.8.b)' 


(Where  I V'/f  I  and  7  were  defined  above  in  connection  with  Eqs.  (3.6)  and  (3.6.a,b,c,d)', 
and  where  cos  4>  =  (3 H/bx')/  W'H  I,  sin  4>  =  (bHjby')l  IV'// 1,  and  the  coordinate  systems 
have  been  defined  such  that  4>  is  equal  to  zero  for  a  perfectly  horizontal  surface.) 


We  have  referred  to  polarization  in  this  section  only  implicitly.  The  subscripts  Aa 
and  Bb  used  in  (3.3 )  through  (3.5)  and  (3. 7.0,..., 4)  refer  to  particular  polarization  situa¬ 
tions.  The  equations  in  Appendix  IV  from  which  (3.7.0,... ,4)  are  justified  contain 
quantities  with  subscripts  VV,  VH,  HV,  and  HH,  where  the  first  subscript  indicates 
received  polarization  and  the  second  subscript  indicates  transmitted  polarization.  Thus, 
e.g.,  if  the  transmitted  and  received  polarizations  corresponding  to  uAa  in  (3.1)  were  both 
vertical,  then  the  subscripts  Aa  would  be  indicated  as  VV.  If  the  transmitted  polarization 
corresponding  to  uBb(t)  were  vertical  and  the  received  polarization  were  horizontal,  then 
subscripts  Bb  would  be  indicated  as  HV,  If  the  transmitted  and  received  polarizations 
were  circular,  then  the  subscripts  Aa  or  Bb  might  indicate  (for  example)  LR,  implying 
that  the  transmitted  polarization  is  right-handed  and  the  received  polarization  is  left-handed. 
In  this  case  (3.1)  would  be  written  as  follows: 


uLR(t)  =  Re\uLR(t)d'JU'ot\, 


(3.9) 


where 


ULR(t)=uVR(t)+jUHR(t), 

where  the  subscripts  VR  and  HR  imply  respectively  vertically  and  horizontally  polarized 
received  radiation  with  right-handed  circular  transmitted  polarization. 


ENSEMBLE  AVERAGES  OVER  SURFACE  FLUCTUATIONS 

The  quantities  used  in  the  final  results  will  be  ensemble-averaged  over  small-scale 
surface  fluctuations.  This  averaging  process  is  discussed  in  Appendix  II  and  results  given 
there  are  used  as  a  basis  for  much  of  what  follows. 
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Averages  of  Voltage  Products 

If  we  expand  uAa(t)  as  given  in  (3.1)  through  fourth  order,  we  obtain 

uAaW  =  UAaoW  +UAal^  +  '  +  uAa4^-  <41> 

The  averaging  of  (3.2)  with  the  aid  of  (3. 7.0,..., 4),  (II.4),  (11.13)  and  (11.16)  leads  to 

~^AaBbo^'  W  +  <“a«B62<*’A*»  +  <4'2> 

where  the  subscripts  0,  2,  and  4  refer  to  perturbation  order  and  where  (it  is  understood  that 
the  argument  of  all  Aa  factors  is  t,  that  of  Bb  factors  is  t  +  At  and  all  quantities  UAaBb 
have  arguments  (f,Af)) 

*  1  -  ju  At 

uAaBBo=~£Re  (UAaBboe  °  )» 

where 

"AaBbo  *  UAaoU*Bbo, 
<'uAaBb2>S‘~2Re^<‘^AaBb2>e  °  )* 

where 

<^AaBb2  >m  UAao<UBb2>*<UAa2>UBbo*<UAalUBbl>' 

(4.2.C)' 

1  jio  At 

<uAaBb4>  ~~^Re  (<UAaBb4>e  °  )• 


(4.2.a)' 


(4.2.b)' 


where 


<^4oBb4>  ~  <UAa2UBb2>  +  <UAal  UBb3>  +  <'UAa2UBb\> 
+  UAao<U*Bb4>  +  <UAa4>U*Bbo- 


(4.2.d)' 
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Averages  of  Power  (or  Quadratically  Rectified  Voltage)  Products 
The  ensemble  averse  of  UAaBb(t,At),  obtained  from  (3.4), 
<UAaRb(tAt)-  =  <UAaAM,0)><0RhRh(t  +  At, 


ble  averse  of  UAaBb(t,At),  obtained  from  (3.4),  is: 
JAaBb  (t,At)-  =  <UAaAa(t,0)XUBbBb(t  +  At,  0)> 

+  <f^aAa^0>-<^a.4a^0)]  - 
^BbBc^  +  ®)“  <'UBbBb(t  +  Af> 


(4.3) 


(*W  <<  +  At •  °>  “  <UBbBb(t  +  Af,  0)>]  >, 

where  each  of  the  factors  in  the  first  term  can  be  considered  as  specializations  of  factors 
in  terms  in  (4.2.a,b,c)'  and  the  second  term,  expanded  to  fourth  order,  with  the  aid  of 
(3.2),  (3.7.0,. ...4),  (11.13)  and  (11.16)  is 

<A^AaB6(^)>  =  <^AaAa^°)-  <^ala«>°»] 


<iuAtAmm-<uAltUm>]m 

BbBb +  ~  BbBb ^  +  At'  ^  = 

^^AaAa  1  (t-0)f/fihRh1  ( t  +  At,  0)> 


^AaAal^^BbBbl^  +  At'  °)> 

+  <'^AaAal^t'^^BbBb3^t  +  At> 

+  <'^AaAa3^t'°)^BbBbl^t  +  At’  °)> 

+  AaAaiS* BbBbZ^  +  At' 

~  <^^AaAa2^t,^'><"^BbBb2^t  +  At'  ®)>* 

whose  individual  terms  are  (where  it  is  understood  that  arguments  of  Aa  facl 
(f,0)  and  those  of  Bb  factors  are  (f  +  At,  0)): 

^AaAal^BbBbl  >-2Re[U*aoUBbo<UAalU*Bbl> 

+  V*AaoU*Bbo<VAalUBbl»> 


(4.4) 


factors  are 


(4.4.a)' 


<U AaAalV BbBbB 


>-2Re[U*aoUBbo<UAalU*Bb3> 


+  UAaoU*Bbo<VAotVBb2> 


+  *-TAt0<'UAalUBblUBb2> 


+  UAao<U*Bbi<alUBb2>h 


(4.4.b)' 
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<l^AaAaZ^BbBbl  >  =  ™e[U*BboUAao<UBblU*Aa3> 


+  UlboV*Aao<UBblVAa3> 


+  U*Bbo<V*Bbl  UAalUAa2> 


+  UBbo<'UAalUBblUAa2>]' 


<UAaAa2UBbBb2 


>  =  2Re[UAnnU*<U*URhi> 


+  UAaoUBbo<UAa2UBb2> 


+  UAao<UBblUBblUAa2> 


+  ^Rba<'^Aa1  ^ Aa\ ^862^1 


+  <V*AalUAalV*BblUBbl> 


'  <UAaAa2><UBbBb2> 


(4.4.c 


(4.4.d 


(4  4 

From  (3.7.0,. ..,4),  (4.2),  (4.2.a,b,c)',  (4.4),  (4.4.a,...,e)\sand  a  number  of  results  in 
Appendi'.  II,  we  can  express  <UAaBb(t,At)>  and  <AUAaBb(t,At)>  in  more  compact 
approximate  forms.  First,  from  (4.2),  (4.2.a,b,c)\  (3.7.0,..., 4),  and  (11.17, 18, 19,  20,  21), 

«^.w(W|V“4,l.  (4.: 

where 

<UAaBb^At)>=UAao<o+UAol<ulb2>^<U*Bb4>] 

+  U*Bbol<UAa2>  +  <UAa4>l 

*  *  (451 
+  <UAaXbl>^<UAa2U*Bb2> 

+  <UAaXb3>  +  <UAaSU*Bbl> 
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where  it  is  understood  that  the  argument  of  UAa  is  always  t  and  that  of  Ugb  is  always 
t  +  At,  and  where  (from  3.7.0) 


27T 


=/ dK  *>  fas  FAa <4-5'2»' 
o  o 

From  (3.7.2)  and  (11.18) 


(4.5.3/ 

(“„■»> 

where  W(k)  is  the  spatial  spectrum  of  the  fluctuations.  From  (3.7.4)  and  (11.21,18,9) 


<^.4>*(— )> 


«V°> 


fSA«4^’  ^l'~hvh2^  +  SAo 4^5’ ^1’ ^2’  -&1  > 

+  SAa4^5>  &1  >  ^2 ’  ~kl M  • 

From  (3.7.1)  and  11.19,14), 

It  271  2 

<uA'A>i>-fafa3  si”«3  fas 

o  o 

SA«l<*o<#3’  ^3’  «  +  AO) 

From  (3.7.2)  and  11.21, 19,14), 

IT 

<UA,A>2>  -  <VA.3><Vlb  2>*  fa)2 fa',  “»«3  - 


(4.5.4)' 


(4.5.5)' 


it 

J *'sfa.K  *3- ‘KiV'r  *'r  ‘  *  wffamkW.. 
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2co 


<-  —  &  .  *3 . 0  -  * *>•  ‘W 


(4.5.6)' 


t^B62^ao^’  ^3’  *  + 


-*  2co 

+  ^B62(^o^3>  ^3>  *  +  ^  2o^3’  ^3>  *  +  ^)  "  £)) » 


where  <UAa2>  and  <UBb2>  axe  specializations  of  (4.5.3)'. 
From  (3.7.1, 3)  and  (11.21, 18,19,14), 

2  7T 


2/ *>3  »in  *J f M'3  K- ~4o«>3- K-‘»- 


FAa^ 3’^3>  f)*W03’^3’  f  + 


(i^MiffdkWiknsl^  (457), 

2c«> 

<&  («$.  ^3*  t  +  At) &  {^,  0' ,  t  +  At)$) 


2(a) 

+  ^B63^o^3’^3’  *  +  Af);fe.  ”^.o^3’^3’  *  +  ^)) 

+  ®B63^o  ^S'^s*  * +  Af);fe*-ft))]. 

From  (3.7.0,. ..,4),  (4.4),  (4.4a, ...e)\  and  (11.17,18,19,20,21,23,24,25,26), 

<AUAaBb(t,At)>^  2  Re  ^laoUBbo<UAayBbl> 

+  u*AaouBho<uAaXv*Bbz>  +  u*BbouAao<vBblu*AaZ> 


(4.6) 


-  ^ao^bo<^a2><^62>-  UAaoU*Bbo<U*Aa2><UBh2> 
*U*nn<URh2><U*n,URh,>^U*<UAai><U*UAa,>), 


where  all  factors  in  terms  of  (4.6)  are  specializations  of  (4.5.2,. ..,7)'. 
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; 


In  the  special  case  where  no  zero-order  terms  exist  (i.e.,  no  specularly  reflecting 
surface  and  no  targets),  then  (4.5)  and  (4.6)  assume  much  simpler  forms.  In  this  case 
(from  (4.5)  with  all  zero-order  terms  set  equal  to  zero) 

<UAaBb^At)>  =  \Re  [<UAalU*Bbl>  +  <UAa2U*Bb2> 

*  (4.7) 

and  (from  (4.6)  with  all  zero-order  terms  set  equal  to  zero) 

(4.8) 


Applications  of  Averages 

Let  us  now  consider  the  applications  of  (4.5).  With  the  aid  of  (4.5.1,.. .7)',  it  can  be 
used  to  obtain  the  following  quantities  of  possible  interest. 

(1)  The  crosscovariance  function  (or  crosscorrelation  function;  abbreviated  as  CCF 
in  what  follows)  between  vertically  (or  horizontally)  polarized  return  signal  voltages  with 
possibly  different  incident  polarizations.  If  the  incident  polarizations  are  the  same  and 
At  =  0,  then  this  quantity  is  the  average  power  in  the  vertically  (or  horizontally)  polarized 
return. 

(2)  The  CCF  between  right-handed  (or  left-handed)  circularly  polarized  signal 
voltages  with  possibly  different  incident  polarizations.  If  the  incident  polarizations  are 
the  same  and  At  =  0,  then  this  quantity  is  the  average  power  in  the  left-handed  (or  right- 
handed)  circularly  polarized  return. 

(3)  The  CCF  between  vertically  and  horizontally  polarized  signal  voltages,  with 
possibly  different  incident  polarizations. 

(4)  The  CCF  between  left  and  right-handed  circularly  polarized  signal  voltages  with 
possibly  different  incident  polarizations. 

Eq.  (4.6)  with  the  aid  of  (4. 5. 2,.. .,7)'  can  be  used  to  determine  the  following 
quantities  of  possible  interest: 

(5)  The  CCF  between  deviations  from  average  power  (or  quadratically  rectified 
voltage)  in  two  vertically  (or  horizontally)  polarized  return  signals  with  possibly  different 
incident  polarizations.  If  the  incident  polarizations  are  the  same  and  At  =  0,  then  this 
quantity  is  the  variance  of  the  vertically  (or  horizontally)  polarized  signal  power. 

(6)  The  CCF  between  deviation  from  average  power  (or  quadratically  rectified 
voltage)  in  vertically  and  horizontally  polarized  returns  with  possibly  different  incident 
polarizations. 
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(7)  The  CCF  between  deviation  from  average  power  or  quadratically  rectified 
voltage  in  two  left-handed  (or  right-handed)  circularly  polarized  returns  with  possibly 
different  incident  polarizations.  If  the  incident  polarizations  are  the  same  and  At  =  0,  this 
quantity  is  the  variance  of  the  left-handed  (or  right-handed)  circularly  polarized  signal 
power. 

(8)  The  CCF  between  deviation  from  average  power  (or  quadratically  rectified 
voltage)  in  left-handed  and  right-handed  circularly  polarized  returns  with  possibly  different 
incident  polarizations. 


Transmit-Receive  Polarization  Modes 

There  are  16  transmit-receive  polarization  modes,  any  one  of  which  may  be  of 
interest  in  a  particular  application.  These  are  given  in  the  chart  below. 


Receive 

polarization 

Transmit 

polarization 

Designation 

(1)  Linear-vertical 

Linear-vertical 

VV 

(2)  Linear-vertical 

Linear-horizontal 

VH 

(3)  Linear-vertical 

Circular-left-handed 

VL 

(4)  Linear-vertical 

Circular-right-handed 

VR 

(5)  Linear-horizontal 

Linear-vertical 

HV 

(6)  Linear-horizontal 

Linear-horizontal 

HH 

(7)  Linear-horizontal 

Circular-left-handed 

HL 

(8)  Linear-horizontal 

Circular-right-handed 

HR 

(9)  Circular-left-handed 

Linear-vertical 

LV 

(10)  Circular-left-handed 

Linear-horizontal 

LH 

(11)  Circular-left-handed 

Circular-left-handed 

LL 

(12)  Circular-left-handed 

Circular-right-handed 

LR 

(13)  Circular-right-handed 

Linear-vertical 

RV 

(14)  Circular-right-handed 

Linear-horizontal 

RH 

(15)  Circular-right-handed 

Circular-left-handed 

RL 

(16)  Circular-right-handed 

Circular-right-handed 

RR 

The  factors  appearing  in  the  terms  of  (4.5)  and  (4.6)  are  all  either  average  voltages 
<I7qg>  or  where  £  =  0,2,  or  4  (the  case  £  =  0,  of  course,  does  not  requ’-e 

ensemble  averaging),  or  averages  of  products  of  two  voltages,  of  the  form  <U j  ,  >  or 

complex  conjugates  of  such  quantities.  We  note  that  in  a  quantity  <UA  g>  (or  *  s), 

A  (or  B)  can  represent  any  one  of  the  16  polarization  cases  listed  in  the  chart.  In  a 
quantity  like  <UA%U%m>,  A  and  B  can  represent  any  of  the  16  cases  listed  in  the  charts; 
hence,  there  are  a  total  of  16  X  16  =  256  possibilities  for  such  a  term.  The  ability  io  choose 
the  desired  case,  of  course,  can  easily  be  programmed  into  a  computer.  The  quantity  SA  g 
appearing  in  (3.7.0,... ,4)  is  calculated  for  Aa  =  VV,  VH,  HV  or  HH.  Extension  of  these 
results  to  include  circular  polarization  can  be  accomplished  through  an  extension  of  (3.9) 
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and  the  observation  that  the  linearly  polarized  response  to  incident  circular  polarization 
is  given  by 


uvrW 

Uyy(t) 

+ 

jUyH 

VL 

vv 

- 

VH 

HR 

HV 

+ 

HH 

HL 

HV 

- 

HH 

The  circularly  polarized  response  to  linearly  polarized  incident  radiation  is 


ULy(t ) 

ii 

c: 

«*e» 

+ 

iUH  V1 

LH 

VH 

+ 

HH 

RV 

VV 

- 

HV 

RH 

VH 

- 

HH 

From  (4.9)  or  (4.10), 

ULL(t)  =  UVL(t)  +]UHL(t)  =  ULV(t)  -  jULH(t) 

=  [Uyy(t)  + 


(19) 


(4.10) 


(4.11.a) 


URrW  =  UVR«)->UHR(t)  =  URV(t)+jURH(t) 

(4.11.b) 

=  l^yy(t)  +  URR(t)]  + 


^lrW  ~  ^vrW  +  Wur^)  ~  ULV(t) + 

(4.11.C) 

/  =  [Uyy(t)~  UHH(t)]  +  j[Uyjj(t)  +  U If  y(t)]  ,  and. 


URL«)  =  UyL(t)-jUHL(t)  =  URV(t)-jURH(t) 

+  [Uyy(t)-  UHH(t )]  -  j[UyH(t)  +  UHy(t)). 


(4.11-d1 


Equations  (4.11.a,...d)  can  be  used  to  form  the  averages  needed  in  (4.5)  and  (4.6)  from 
the  basic  calculated  quantities  Uvv,  UVH,  URV,  Ujjfj. 


RAEMER 
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Inside  the  integral  in  each  of  these  quantities  (see  3. 7.0,. ..,4)  is  the  quantity  SAa%, 
where  Aa  =  W,  VH,  HV  or  HH.  On  the  computer,  cases  such  as  LL,  LR ,  etc.  will  be 
implemented  by  programming  the  statement 


SLLi  ~  (SVV9.  +  SHHi^  +  )(SHVi  ~  SVHt^‘ 
RR  fi  + 

LRi  +  + 

RLi  + 


(4.12) 


This  is  the  easiest  way  to  handle  cases  involving  circular  polarization  on  both  trans¬ 
mission  and  reception.  The  averages  required  in  (4.5)  and  (4.6)  could  be  (and  have  been) 
calculated  explicitly  for  these  cases  but  most  of  the  expressions  obtained  for  them  arc  very 
long,  are  net  needed  in  subsequent  calculations,  and  do  not  lend  themselves  to  easy  physical 
interpretation.  Hence,  they  will  not  be  included  in  this  report. 


RESULTS  WITH  FIRST-ORDER  FIELDS  ONLY 

In  this  section  we  will  summarize  and  discuss  the  results  obtained  in  Appendix  IV 
for  the  case  where  only  the  first-order  term  in  the  polarization  matrix  is  considered. 
Specular  reflection  terms,  target  returns  and  second,  third  and  fourth  order  terms  are  not 
included. 

We  can  obtain  a  certain  amount  of  information  from  these  results  but  unfortunately 
cannot  really  address  the  problem  of  cross-polarized  components.  However,  a  certain 
“quasi-cross-polarization”  effect  exists  and  can  be  studied  through  these  results.  This 
point  will  be  elaborated  upon  later. 


Polarization  Matrix  Element 

The  polarization  matrix  element  corresponding  to  a  [1]  is  calculated  in  Appendix  IV. 
The  result,  given  in  Appendix  IV,  is 


VH 

KV 

HH 


v  v 

H  H 
V  H 


-  ZlozKz+X 


ozHoxZ!v>z(i 


->0 


V 

H 

H 


t 


(5.1) 
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(Note  that  the  first  subscript  V  or  H  refers  to  received  signal  polarization  and  the  second 
subscript  refers  to  transmitted  signal  polarization,  where  V  and  H  denote  vertical  and 


horizontal  respectively.) 

The  quantities  in  (5.1)  are  defined  below: 

v  =  complex  refractive  index  of  the  medium  (5.2.a) 

oo  =  (angular)  radar  frequency  (5.2.b) 

c  *  free  space  velocity  of  light  (5.2.e) 

A<  bo!)  =  < y0z  -  “o*)*1  -  alz  -  aoz Toz>  <5-2d> 


9io  =4r’  = 


(from  I.27.a)  (5.2.f) 

olqx  =  sin  0  g  cos  0g  (cos  6  cos  7  cos  $  -  sin  6  sin  y) 

+  sin  0 3  sin  0g  (-  cos  5  sin  4») 

+  cos  0g  (cos  8  sin  7  cos  $  +  sin  8  cos  7);  (from  1.28) 
otoy  =  sin  0g  cos  0g  (-  cos  7  sin  4>)  +  sin  0g  sin  0 g  (-  cos  d>) 

+  cos  0g  (-  sin  7  Rin  4>):  (from  1.28) 

■  W 1,2  -  +aoy) 

^oz 

where  aox  and  aoy  are  defined  in  (5.2.g.x,y) 

7  7  angle  of  peak  of  antenna  beam  relative  to  the  vertical  (5.2.h) 


(5.2.g.x) 

(5.2.g.y) 

(5.2.g.z) 


a  = 
-*o 


a 


ox 


a 


oy 


a 


ox 


oy 


a 


oz 


(from  1.22) 


(5.2.e) 
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(^3^3)  -  spherical  angles  of  the  illuminated  point  on  the  surface  relative  to  the  radar  in 
the  S'3  system  (whose  z'3  l:  is  is  along  the  peak  of  the  antenna  beam) 


.  J. 

cos  0  -  ) - - p-  (from  III.15.a) 

(5.2.j.l) 

„  iv'/f  1 

sm  6  =  , - 1 - r  (from  I11.15.b) 

V^l  +  )v'H  j2 

(5.2.j.2) 

H  ' 

cos  4>  -  (from  III.lS.c) 

IV  H 1 

(5.2.J.3) 

Hy' 

sin  4>  -  (from  III.15.d) 

IV  H 1 

(5.2.j.4) 

where 

H  -  H(x\  y)  =  Mean  surface  as  a  function  of  (x1,  y'),  the  horizontal 
coordinates  in  the  S'  system  (whose  z  axis  is 
perfectly  vertical  and  whose  (x\  y')  plane  is  along 
the  horizontal  sea  surface) 

(5.2.j.5) 

|)CJ 

»,•  - 

(5.2.j.6) 

y  dy  V  ,y  ’ 

(5.2.j.7) 

V'H  =v4ffx')2  *  (H/)2 

(5.2.j.8) 

=  2v*  "  ~£V*ao*  (frora  128  wthj  =  -&o,  pz  =  -ctoz) 

Vy  Vy  oz  V  oy 

Hx  Hx  H  ox 

Hy  Hy  H  oy 

(5.2.k) 

=  cos  8  cos  7  cos  $  -  sin  5  sin  7 

(5.2.S.1) 

£Vy  =-cos7  6in<?> 

(5. 2. £.2) 

=  sin  5  cos  7  cos  <I>  +  cos  8  sin  y 

(5. 2. £.3) 
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lHx  =  -  cos  6  sin  4> 

(5.2.C.4) 

y  =  -  cos  4> 

(5.2.C.5) 

iHz  =  -  sin  5  sin  $ 

(5.2.C.6) 

(Eqs.  (5.2.£.1,...,6)  are  from  (1.27)) 

The  vector  and  scalar  products  in,  (5.1)  are  as  follows: 


\  =  “  =  ^Hx^Vy  ~  ^Hy^Vx 

(5.3.a) 

=  (£hxShK  =  0 

(5.3.b) 

(%oX$y) 2  =  “oxCVy  "  “oyVx 

H  H  H 

(5.3.c) 

OL  *2 j/  =  ct  £»»  +  otni'v 
-*o  -►v  ox  Vx  oy  Vy 

H  H  H 

(5.3.d) 

In  the  degenerate  case  where  the  mean  surface  is  perfectly  horizontal,  some  of  the 
quantities  defined  in  (5.2.a,...C)  and  (5.3.a,...,d)  are  significantly  simplified.  In  this  case 

H(x',y')  =  0. 

(5.4) 

Choosing  to  set  Hy<  identically  to  zero  and  allowing  Hx,  to  become  arbitrarily  small 
(just  an  arbitrary  choice  of  coordinate  orientation  which  does  not  affect  the  end  results),  we 
obtain  from  (5.2.j.l,... ,4): 

cos  4>  =  1 

(5.5.a) 

sin  $  =  0 

(5.5.b) 

IV'tfl  =  0 

(5.5.C) 

cos  8  =  1 

(5.5.d) 

sin  S  =  0 

(5.5.e) 

Using  (5.5.a,...,e)  in  (5.2.g.x.y.z),  (5.2.k),  and  (5.2.£.6),  we  obtain 

olQx  ~  sin  0g  cos  <f> g  cos  y  +  cos  0g  sin  y 

(5.6.a) 

aoy  =-sin  0g  sin  0g 

(5.6.b) 
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oz 


y 


oz 


=  + 


1 

-  [sin2  0g  (cos2  0g  cos2  7  +  sin2  <f>3)  +  cos2  8 g  sin2  7  ... 

,2 

(5.6.c) 

+  2  sin  0g  cos  d'3  cos  <t>'3  cos  7  sin  7]  l 


lVx  =  cos  7 

(5.7.a) 

*=> 

11 

0 

(5.7  .b) 

iVz  =  sin  7 

(5.7.c) 

0 

11 

X 

0? 

(5.7.d) 

=  “  1 

(5.7 .e) 

=  0 

(5.7  .f) 

To  check  our  results  against  those  previously  obtained  by  other  workers,  we  further 
specialize  to  the  case  where  the  surface  is  perfectly  horizontal  and  illumination  is  in  the 
(X3  -  Z3 )  plane  and  along  the  peak  of  the  beam,  corresponding  to  the  conditions 

03=0 


Using  (5.8)  in  (5.6.a,b,c),  (5.2,k),  (5.7.a,...,f)  and  (5.3.a,b,c,d),  we  obtain: 


aox =  sin  y 

(5.9.a) 

0 

n 

(5.9.b) 

=  -  y/l  -  sin2  7  =  -  cos  7 

(5.9.c) 

yQz  - sin2  y 

(5.9.d) 

iVx  =  cos  7 

(5.9.e) 

11 

0 

(5.9.f) 
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iVz  =  sin  7 

CWx  =  0 
eHy  =  ~  1 
V-Hz  =  0 
i'Vx  =  sec  y 

fi'vy  =  0 
*//*  =  0 

(%X^Z  ~  0 

(«o‘^)  =  tan7 

(&>*&>* 

(£>*&> =  0 

(&*^)2  =  sec  7 

A(  Mo  l> =  (T „*  "  %*)<!  "  «022  “  «0*  7oa)  =  I\A2  -  sin2  7  +  cos  7] 
[sin2  7  +  cos  y\Jv 2  -  sin2  7]  =  p2  cos  7  +y/v2  -  sin2  7. 
Substitution  of  (5.9.a,...,t)  and  (5.3.b)  into  (5.1)  results  in: 

2/^(p2  ~  1)  cos2  7[(p2  -  1)  sin2  7  +  v 2] 


ISvy(~  ^0)34J(i ) = 


(p2  cos  7  +y4>2  -  sin2  7)2 


(5.9.g) 
(5.9.h) 
(5.9.i) 
(5.9.j) 
(5.9.k) 
(5.9.8) 
(5.9  .m) 
(5.9  .n) 
(5.9.0) 
(5.9.p) 
(5.9.q) 
(5.9.r) 
(5.9.s) 

(5.9.t) 

(5.10.a) 

(5.10.b) 


(5.10.C) 
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“2 j£(v2  -  1)  COS  7 

t SHH^~  ao^a(1)=  rz . ,  2  Z,  •  (5.10.d) 

(cos  7  +y/v*  -  sinz  yy 

The  results  given  in  (5.10.a,b,c,d)  are  consistent  with  those  reported  in  the  literature, 
e.g.  Barrick  and  Peak  [46] ;  Guinard  and  Daley  [10] ;  Guinard,  Ransone  and  Daley  [14] ; 
and  Valenzuela,  Laing  and  Daley  [13] .  Notation  and  definitions  of  angles  differ  from  those 
used  here  but  when  the  meanings  of  the  symbols  in  those  papers  were  interpreted  in  terms 
of  the  author’s  notation  and  definitions,  the  results  were  found  to  be  in  exact  agreement 
with  those  reported  previously. 

Let  us  refer  back  to  the  definition  of  vertical  and  horizontal  polarization  used  in  this 
report.  We  define  vertically  polarized  transmission  or  return  (labeled  V)  as  that  whose  E 
vector  is  in  a  given  direction  on  the  antenna  aperture  and  horizontally  polarized  transmission 
or  return  (labeled  H)  as  that  whose  E  vector  is  perpendicular  to  the  V  direction.  But  the 
E  vector  corresponding  to  a  plane  wave  of  either  V  or  H  type  (according  to  the  above 
definition)  launched  at  the  aperture  and  striking  the  sea  surface  at  a  given  point  has 
components  both  vertically  and  horizontally  polarized  (according  to  the  standard  defini¬ 
tion  of  vertically  polarized  radiation  as  that  whose  vector  is  in  the  plane  of  incidence 
and  horizontally  polarized  radiation  as  that  whose  JE£  vector  is  normal  to  the  plane  of 
incidence)  and  these  components  are  acted  upon  differently  by  the  surface  reflection. 

Hence  the  return  corresponding  to  a  V-type  transmitted  wave  has  an  H-component  and  the 
return  corresponding  to  an  H-type  transmitted  wave  has  a  V-component.  If  the  return  is 
considered  as  a  superposition  of  plane  waves  from  all  angles  contained  within  the  antenna 
beam  pattern,  then  each  such  plane  wave  has  both  V  and  H  components  regardless  of 
whether  their  corresponding  transmitted  waves  were  of  V  or  H  type.  The  total  2£  field  on 
the  aperture  resulting  from  this  superposition  of  returns  has  a  different  ratio  of 
V-component  to  H-component  than  that  resulting  from  the  superposition  of  ^-fields  of  all 
the  transmitted  waves.  This  constitutes  a  kind  of  “quasi-cross  polarization”  effect  which 
appears  in  the  return  even  in  first  order,  although  it  is  very  well  known  that  no  real  cross 
polarization  (defined  in  the  standard  way,  i.e.,  referred  to  polarization  at  the  sea  surface, 
not  at  the  antenna  aperture)  exists  in  first  order  return. 

The  exception  to  this  is  the  case  where  all  of  the  radiation  goes  out  in  a  single 
direction,  i.e.,  along  the  antenna  beam  axis,  in  which  case  the  components  of  the  ^-vector 
of  the  transmitted  or  returned  plane  wave  correspond  to  the  components  of  the  ^-vector 
on  the  antenna  aperture.  This  is  the  case  covered  by  Eqs.  (5.8),  (5.9.a,...,t)  and  (5.10.a,b,c,d). 

As  soon  as  a  small  departure  from  this  simple  situation  occurs,  then  the  return  has  the 
feature  indicated  above,  i.e.  its  “polarization”  as  defined  at  the  antenna  aperture  is 
different  from  that  of  its  transmitted  wave. 

This  is  true  even  for  a  perfectly  horizontal  mean  sea  surface.  Of  course,  the  wider  the 
beam,  the  greater  is  this  “quasi-cross  polarization”  effect,  because  a  wider  beam  has  much 
more  wave  energy  propagated  at  angles  far  off  the  beam  axis.  With  a  narrow  beam,  e.g. 
only  a  degree  or  two  wide,  there  should  be  very  little  of  the  effect,  because  at  small 
angles  off  of  the  beam  axis,  the  J£-field  components  defined  as  vertically  and  horizontally 
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polarized  at  the  antenna  aperture  come  much  closer  to  being  the  vertically  and  horizontally 
polarized  components  at  the  reflecting  surface. 

For  any  given  beamwidth,  “quasi-cross  polarization”  effect  would  be  further  enhanced 
by  a  rough  mean-surface,  because  waves  propagating  out  at  different  angles  would  encounter 
a  wide  range  of  angles  of  incidence  at  the  surface  and  hence  the  probability  of  a  vast 
difference  between  the  polarization  defined  at  the  antenna  and  that  defined  at  the  surface 
would  be  much  greater,  even  at  angles  not  far  off  the  beam  axis. 

To  summarize  this  point,  our  concept  of  polarization,  which  defines  vertical  and 
horizontal  polarization  in  terms  of  two  mutually  perpendicular  directions  of  the  f£-field 
on  the  antenna  aperture,  does  not  necessarily  correspond  to  the  standard  definition  of 
vertical  and  horizontal  polarization,  which  relates  to  the  §-field  in  the  plane  of  incidence 
or  normal  to  the  plane  of  incidence  respectively.  Suppose  we  think  of  the  return  as  a 
superposition  of  plane  waves  reflected  from  the  random  surface.  If  the  mean  surface  is 
horizontal,  the  only  one  of  these  waves  for  which  these  two  definitions  coincide  is  the  wave 
along  the  beam  axis.  In  the  case  in  which  the  results  (5.10.a,b,c,d)  are  obtained,  that  is  the 
only  wave  propagating;  cross  polarization  terms  vanish.  It  can  easily  be  shown  that  the  VH 
and  HV  terms  also  vanish  (with  a  perfectly  horizontal  mean  surface)  in  the  case  where 
the  illumination  is  in  the  (*3  -  23)  plane  (0 3  =  0).  In  this  case,  a  V  type  wave  ( our  defini¬ 
tion)  has  no  components  that  are  horizontally  polarized  ( standard  definition)  at  the  sea 
surface  and  an  H-type  wave  (our  definition)  has  no  components  vertically  polarized 
(standard  definition)  at  the  sea  surface.  This  is  a  degenerate  case  in  our  analysis  because  it 
covers  only  one  plane  of  illumination  and  we  wish  this  analysis  to  cover  a  superposition 
of  waves  from  an  arbitrary  solid  angle. 

To  show  the  vanishing  of  VH  and  HV  terms  in  this  case,  we  need  only  observe  that 
setting  03  equal  to  zero  in  (5.6.a,b),  and  (5.7.a,b,d,e)  and  invoking  (5.3.b)  results  in: 


«0*  =  sin  (0g  -  7) 


(5.11.a) 


a  =0 
oy 


(5.11.b) 


i'vx  =  cos  7  - 


sin  7  sin  (0g  -  7) 


(5.11.C) 


(5.11.d) 


(5.11.e) 


=  "  1 


(S.ll.f ) 


(v!vMi1'v)=0 


(5.11.g) 


1  i £*Sfc  « 
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<511h) 

Eqs.  (5.3.b)  and  (5.11.g,h)  imply  that 

ISvHl-LKai  D=  fS//v(-io)Ja(D  =  °  <512) 

for  the  special  case  where  0g  =  0. 


Averages  Involving  Only  First-Order  Terms 

The  generic  average  of  the  product  of  two  polarization  matrix  elements  each  of  which 
contains  only  first-order  terms  is  the  term  <UAal  U%b\>  “  (4.2.b)\  which  is  of  second 
order.  The  generic  average  of  power  (or  quadratically  rectified  voltage)  products  in  the 
case  where  the  polarization  matrix  contains  only  first-order  terms  is  the  term 
^^Aal^Aal^Bbl^Bbl^  'n  (4.4.d)',  whi«*h  is  of  fourth  order. 

Summarizing  the  applicable  results  in  the  previous  section,  wc  have  (from  (4.2.b)') 

=  <«•«> 

and  from  (4.8) 

<WAaBb (t,Af)>-  |<tfAal^61>K  (5.14) 

where  (from  (4.5.5)'  with  a  slight  notational  change) 

7T  2j r 

<'UAa\UBbl>  =  (~^}2 J  d03  sin  J  d<l> 3  ^(“  ~  %  W 3»03»O)- 

0  0  (5.15) 

FAM*3*)*ZbVa***  +  W 

lSAa(-  «o(03^3>f))]a(l)[SB6(-i,(034,f  +  Af))]0*(l). 

Here  FAa  and  FBb  are  given  by  (3.5);  the  various  constituents  of  FAa  and  FBb  are  in 
tum  obtained  from  (3.6)  and  (3.6.a,b,c,d)'.  The  quantities  [SAa]  and  [SBb]  are  of  course 
the  polarization  matrix  elements  (5.1),  or  combinations  thereof.  The  pairs  of  subscripts 
(Aa)  or  (Bb)  on  the  [S]  's  may  be  (VV),  (VH),  (HV)  or  (HH)  if  attention  is  confined  to 
linearly  polarized  transmission  and  reception.  They  may  also  denote  (LL),  (LR),  ( RL )  or 
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The  subscripts  (/la)  and  ( Bb )  can  each  denote  any  of  the  possibilities  covered  by 
(5.16),  (5.17)  and  (5.18)  as  well  as  the  standard  cases  involving  only  linear  polarization 
in  both  transmission  and  reception. 


LISTING  OF  SELECTED  RESULTS 
(in  order  of  appearance) 

1..  General  matrix  expressions  for  fields  of  zero-order;  Eqs.  (2.16.0),  (2.16.0)',  p.  19. 

2.  General  matrix  expression  for  fields  of  first  order;  Eq.  (2.16.1),  p.  19;  second  order; 
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Appendix  I 

COORDINATE  SYSTEMS 

The  basic  coordinate  system  used  in  solving  the  boundary  value  problem  is  a  right- 
handed  rectangular  system  S  =  (x,y,z),  with  basis  vectors  (  i^x .  i*y,  i^).  The  system’s  (x,y) 
plane  is  tangent  to  the  tilted  mean  surface  at  its  origin  of  coordinates,  denoted  by  0.  We 
deiine  another  coordinate  system  S'  =  (x',y',z')  with  basis  vectors  (]^x',  i^z'),  its 

(x\  y')  plane  along  the  horizontal  “perfectly  calm”  sea  surface,  its  origin  at  an  arbitrary 
point  0'  on  that  surface,  and  its  z  axis  vertically  upward.  We  are  interested  in  the  trans¬ 
formation  between  these  two  systems. 

Given  that  the  origin  0  of  the  (*,y,z)  system  is  located  at  a  point  (x0',  y0',  z0')  in  the 
(x\  y',  z')  system  (see  Fig.  1-1),  we  can  construct  an  intermediate  system  S"  =  ( x",  y",  z") 
with  basis  vectors  (Xx'>  Xy"<  Xz'X  whose  coordinates  are  parallel  to  those  of  system  S' 
and  whose  origin  is  at  0.  A  point  with  coordinates  (*',  y',  z)  in  system  S'  will  have 
coordinates  (*'  -  x0',  y'  -  y0',  z  -  z0')  in  system  S",  i.e., 


O’ 


Fig.  1-1  —  Systems  S,  S',  s" 
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Referring  to  Fig.  1-2,  we  construct  a  second  intermediate  system  S' "  =  (x' ",  y' ",  z  ") 
with  origin  at  0,  whose  z  "  axis  is  parallel  to  the  z"  axis  of  system  S"  and  that  is  rotated 
by  an  angle  <l>  with  respect  to  S"  such  that  the  x' "  axis  is  in  the  direction  of  the  iocal 
tilt,  i.e.,  by  rotation  in  the  (*' ",  z  ")  plane  we  will  be  able  to  transform  to  the  system  S. 

The  transformation  between  S' "  and  S"  is  (from  Fig.  1-2) 

x'"  =  x"  cos  <1>  +  y”  sin  4>, 

y'"  =  -  x”  sin  4>  +  y"  cos  <l>,  and  (1.2) 


To  transform  from  S'  "  to  S,  we  now  tilt  along  the  (x' " .  z'  ”)  plane  (see  Fig.  1-3) 
through  an  angle  5,  leading  to 

x  -  x'"  cos  8  -  z' "  sin  6, 

y  =  y'",and  (1.3) 

_  t  tt  »  C  .  Iff  c 

z  =  x  sin  6  +2  cos  5. 

: 

i 

i 


NRL  REPORT  8369 


To  obtain  the  transformation  between  S'  and  S,  we  first  use  (1.2)  and  (1.3)  to  express 
(x,y,z)  in  terms  of  ( x ",  y",  z")  resulting  in 

x  =  x "  cos  5  cos  4>  +  y"  cos  8  sin  «l>  -  z"  sin  8, 

y  =  -  x"  sin  <!’  +  y"  cos  <I>,  and  (1.4) 

z  =  x"  sin  6  cos  <I>  +  y"  sin  8  sin  *i»  +  z"  cos  8. 

We  then  use  (1.1)  and  (1.4)  to  obtain  (x,  y,  z)  in  terms  of  (x',  y\ z),  as  follows: 

x  =  (x  -  xc')  cos  8  cos  <1>  +  (y'  -  yQ')  cos  8  sin  <I>  -  ( z '  -  zQ')  sin  8, 

y  =-(*'-  xQ')  sin  4>  +  (y'  -  yQ')  cos  <l>,  and  (1.5) 

z  -  (x'  -  xQ')  sin  6  cos  <1>  +  (y'  -  yQ')  sin  8  sin  <!>  +  ( z '  -  zQ ’)  cos  8. 

Inversion  of  (1.5)  is  straightforward  and  leads  to 

x'  =  Xq  +  x  cos  8  cos  -  y  sin  <!>  +  z  sin  8  cos  4>, 

y‘  =  '  +  x  cos  8  sin  <!>+>'  cos  4>  +  z  sin  8  sin  4>,  and  (1.6) 

z'  -  zQ’  -  x  sin  8  +z  cos  8.. 

It  is  important  to  have  the  transformation  between  the  sets  of  basis  vectors  ( 4* .  Xy  > 
X*)  and  (X*\  Xy’>  Xx’)-  These  are  obtained  by  writing 


x  =Sxx  +iyy  +J*z  -  xo>  +iy ,(yl  -  yo>  +iz(z'  ~  2o> 


(1.7) 
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aid  then  using  (1.5)  or  (1.6)  to  solve  for  (_f  x,  \y,  m  terms  of  (i^x\  ±y',  _£,'),  or  vice 
versa.  The  results  are 

i  =  i  '  cos  5  cos  <!>  +  *'/  cos  6  sin  4>  -  i  1  sin  5, 

~*x  -*x  -►y  -*z  y 

»‘  =  -  i„ '  sin  <l>  + »'  1  cos  4>,  and  (1.8) 

—+y  -+y 

i ,  =  /  '  sin  8  cos  4>  +  i  '  sin  5  sin  4>  +  i  '  cos  5, 

-+z  -+X  ->.y  ->z  ’ 

or  the  inverse  of  (1.8), 

i  J  =  i*  cos  6  cos  4>  -  i  sin  4>  +  i  sin  6  cos 4>, 

-►x  -*x  -*y  -+2 

A  '  =  ir  cos  6  sin  <l>  +  i  cos  4>  +  /  sin  5  sin  <t>  and  (1.9) 

— ►  y  —tx  -►/  -►  z 

i  , '  =  -  sin  5  +  i  cos  8 

We  now  construct  another  coordinate  system  S2'  =  (x2\ y2  ,  z2)>  whose  (x2\  y2) 
plane  is  parallel  to  the  (*',  y')  plane  of  system  S',  whose  origin  is  at  the  radar,  and  whose 
z2  axis  is  vertically  downward.  (See  Fig.  1.4.) 


-6- 

(Ty- 


V2‘  O 


Radar  Position 


<xr’.  Yr'.  *r‘> 


X* 


Fig.  1-4  —  Systems  s',  S^' 


The  transformation  between  S2'  and  S'  is 

Xg  —  —  Xr  +  X  , 


y2  =  yR'-y',  and 


z2  ~  ZR  “ 2  i 

where  ( xR ',  yR ',  zR ')  are  the  radar  coordinates  in  S'. 


(UO) 
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Fig.  1-5  —  Systems  S2\  S 3' 


or 


x2  =  *3’ cos  T +z3>  sin  7* 


y2'=  y3>md 

z2  ~  ~  *3'  sin  7  +  z3'  cos  7 


and  for  the  unit  basis  vectors 


lx  =  lx  cos  7  "  K  sin  7, 

x3  x2  Zz 


iv  ,  =  iv  '.and 


■*y  3  ~"y2 


^3  =^2  SinT  +  S  COS7’ 


or 


i*a  =i*3  cos7+i,23  sin 7, 


VW’,nd 

i  '  =*  -  1'  '  sin  7  +  i  '  cos 


/ 
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From  (I.16.a,b),  (1.14),  and  (1.8)  or  (1.9)  we  obtain 


i^x  =  i>x  '  (cos  5  cos  7  cos  -  sin  5  sin  7)  +_^v  '  (-  cos  5  sin  4>) 


+  i  '  (cos  6  sin  7  cos  4>  +  sin  5  cos  7) 


i^y  =  i^x  '  (-  cos  7  sin  <l>)  +^v  '(-cos4>)+^?  '  (- sin  7  sin  4>)>  and  (I.17.a) 


i  =  i  '  (sin  5  cos  7  cos  <J>  +  cos  5  sin  7)  +  1  '  (-  sin  5  sin  <I>) 


+  1  '  (sin  5  sin  7  cos  -  cos  6  cos  7), 

-*•3 

or 

i  =i  (cos  6  cos  7  cos  -  sin  5  sin  7)  +  i  (-  cos  7  sin  4>) 

->  X  g  — >■  -*• 

+  1  z  (sin  5  cos  7  cos  +  cos  b  sin  7), 

/  '=/'(-  cos  6  sin  4>)  +  1  (-  cos  <I>)  +  1  (-  sin  5  sin  <t>),  and  (I.17.b) 

— >•  y  ^  — >  y  — > 

i ,  '  =  i„  (cos  §  sin  7  cos  +  sin  6  cos  7)  +  i  (-  sin  7  sin  4>) 

— yz  ^ 

+X  z  ^s*n  ^  s*n  7  cos  ^  “  cos  ^  cos 

From  (I.17.a,b)  we  can  obtain  the  relationship  between  the  (x,y,z)  components  and 
(x3',  y3',  z3')  components  of  an  arbitrary  vector.  We  begin  with  the  observation  that 


=  a  i  +  a  i  +  a  i  =  a  '  i  '  +  a  '  i  '  +  a  '  i 

X  -*x  y-*y  Z-+Z  *3  -*xg  y3  ->y  3  z3  -+z  3 

Substituting  (I.17.a)  into  the  LHS  of  (1.18),  we  obtain  the  matrix  equation 


(1.18) 


$ 

i 

i 

1 

r  “i 

t 

ax 

(cos  5  cos  7  cos  $  (-  cos  7  sin  4>)  (sin  5  cos  7  cos  ‘F 

X3 

-  sin  5  sin  7)  +  cos  5  sin  7) 

t 

°y 

S3 

= 

(-  cos  5  sin  <l>)  (-  cos  4>)  (-  sin  5  sin  4>) 

f 

(cos  8  sin  7  cos  $  (-  sin  7  sin  <h)  (sin  5  sin  7  cos  4> 

1 

3 

+  sin  6  cos  7)  -  cos  5  cos  7) 

(I.19.a) 
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Substituting  (I.17.b)  into  the  RHS  of  (1.18)  or  equivalently  inverting  the  3  X  3  matrix 
in  (I.19.a),  we  obtain 


(cos  6  cos  7  cos  'J*  (-  cos  5  sin  4>)  (cos  8  sin  7  cos  <1> 

-  sin  5  sin  7) 


(-  cos  7  sin  <b)  (-  cos  4>) 


+  sin  6  cos  7) 
(-  sin  7  sin  4>) 


(sin  5  cos  7  cos  4>  (-  sin  8  sin  <J>)  (sin  5  sin  7  cos  4> 


+  cos  5  sin  7) 


cos  5  cos  7) 


k; 


,(I.19.b) 


It  is  of  interest  to  have  the  transformations  between  the  (x,y,z)  system  and  the 
(*3',  y3',  23')  system.  To  develop  these  transformations,  we  first  invoke  (1.12)  and 
(I.15.b)  to  obtain 


X 

(cos  7  cos  5  cos  <t>  (-  cos  8  sin  4>)  (sin  7  cos  6  cos  4> 

V 

-sin  7  sin  5)  +  cos  7  sin  6) 

y 

= 

(-  cos  7  sin  4>)  (-  cos  <f>)  (-  sin  7  sin  4») 

V 

z 

(cos  7  sin  8  cos  (-  sin  5  sin  $)  (sin  7  sin  6  cos  <I> 

2  3 

«J 

+  sin  7  cos  6)  -  cos  7  cos  5) 

J  (1.20  .a) 


(cos  8  cos  <I>)  (cos  8  sin  4>)  (-  sin  6) 

(-sin4>)  (cos4>)  0 

(sin  8  cos  <t>)  (sin  5  sin  <I>)  (cos  5) 
Then  from  (I.15.a)  and  1.13),  we  obtain 


( xr  *0  ) 

(Xr'-V) 
(zr  _z0  ) 


V 

(cos  7  cos  5  cos  4>  (-  cos  7  sin  (cos  7  sin  5  cos  4> 

X 

-sin  7  sin  6)  +  sin  7  cos  5) 

y3' 

= 

(-  cos  8  sin  4>)  (-  cos  <t>)  (-  sin  5  sin  4>) 

y 

(sin  5  cos  5  cos  4>  (-  sin  7  sin  4>)  (sin  7  sin  5  cos 

z 

+  cos  7  sin  8)  -  cos  7  cos  5) 

-  cos  7 
0 

-  sin  7 


0 

1 

0 


-  sin  7 
0 

cos  7 


—  — 

(XR  -  Xq  ) 

(yR  ~  y0  ) 

(ZR  ~  z0  ) 
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The  unit  basis  vectors  along  the  spherical  coordinate  axes  in  the  S3'  system  are 
given  by 

i  '  =  i  '  sin  0o '  cos  <£  '  +  i  ,  '  sin  OJ  sin  <f>J  +  i  'cos  Or.', 

~*x  3  3  3  3  3  3  Z3  3 

in  '  =  i  '  cos  0^  cos  <t>J  +  i  '  cos  0, '  sin  0, '  -  i  'sinO„',  and  (1.21) 

-*■”  3  ^3  d  d  ~*y  3  3  3  -**3  3 

i  .  '  =  -  i  '  sin  4>J  +  i  '  cos  <p.J. 

-*v 3  -*X  3  6  -*y  3  3 

Substituting  (I.17.b)  into  (1.21)  we  have 

i  '  -  i  [sin  6^'  cos  <j>„'  (cos  5  cos  7  cos  <1>  -  sin  5  sin  7) 

->-r3  -*x  '■3  6 

+  sin  d3'  sin  03'  (-  cos  6  sin  <1>)  +  cos  03'  (cos  5  sin  7  cos  ‘1> 

+  sin  5  cos  7)]  +  _£  [sin  03'  cos  03'(-  cos  7  sin  <1>) 

+  sin  03'  sin  <t>3'(~  cos  4>)  +  cos  03'(-  sin  7  sin  d’)]  (1-22) 

+J^[sin  03'  cos  <p3'(sm  5  cos  7  cos  <I>  +  cos  5  sin  7) 

+  sin  03'  sin  <P3'(~  sin  5  sin  <1>) 

+  cos  03'(sin  5  sin  7  cos  4>  -  cos  5  cos  7)] . 

Defining x3'  and  y3'  as  the  electric  field  directions  for  vertically  and  horizontally 
polarized  waves  respectively,  and  denoting  x3'  and  y3'  components  with  subscripts 
V'  and  H  respectively,  we  can  write 

EjV%u) = yT  $(§,<•>)  =y  ■  I.  - and 


(1.23) 


EiH{^)  =  iHTEi(B,u)^H’El 


(where  the  scalar  products  on  the  left  are  in  vector-matrix  notation  and  those  on  the  right 
are  in  standard  vector  notation)  where  (from  (I.19.a)) 


(cos  6  cos  7  cos  $  -  sin  5  sin  7) 


Cy  =  evy  =  (-  cos  7  sin  <f>) 


(sin  5  cos  7  cos  4>  +  cos  5  sin  7) 


(I.24.a) 
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x  (-  cos  8  sin  4>) 

%H  =  8Wy  =  (-cos<I>) 

£Wa  (-  sin  5sin  <1>) 

Using  (2.8)  in  (1.23),  we  obtain  the  matrix  equation 


«V,£rV,)  Eu 


«».  (V  ±rtH;»  £, 


(1.24. b) 


(1.25) 


where  and  £Wjc  are  given  in  (I.24.a,b). 


Inversion  of  (1.25)  results  in 


Eiy  I^Vjc  *^3  ^Vz^Hy  ±p  &Hz^~  ^Hx  ±q  ^Hz^Vy  ±p  CVz^ 


Py  En 

^Hy  ±q  * Hz )  “^Vy  ±q  8Vz*  r 


(1.26) 


i  ; 

i  1  '• 


'^Hx  ±o  *Hz)  ±a  C Vz )  ,W 

r 


( 
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The  matrix  H  introduced  in  (2.18.a,b)  is  (from  (1.24.a,b)) 


1 

Hx 

(cos  5  cos  7  cos  <l> 

-  sin  5  sin  7) 

(-  cos  6  sin  <I> ) 

1= 

h'y 

Hy 

= 

(-  cos  7  sin  (I>) 

(-  cos  'I’) 

1 - 

<: 

2 Hz 

(sin  5  cos  7  cos  <l» 

+  cos  6  sin  7) 

(-  sin  6  sin  <l>) 

(1.27) 


The  matrix  L(c£0)  appearing  in  (2.31.  0,1,..„4)  is  obtained  from  L((J) as  given  in  (2.19) 
which  is  equivalent  to  the  2X2  matrix  in  (1.26),  where  =  aQ  and  is  obtained  from 
i  '  (given  by  (1.22)),  i.e.. 


a  = 
-*  o 


a 


a 


oy 


(I.27.a) 


where 


aQX  =  sin  0gs'  cos  03S'  (cos  6  cos  7  cos  «l>  -  sin  6  sin  7) 

+  sin  0 3S'  sin  03S'  (-  cos  5  sin  «l>) 

+  cos  (?3S'  (cos  5  sin  7  cos  ‘h  +  sin  6  cos  7),  and  (1.28) 

uoy  =  sin  03S'  cos  03S'  (-  cos  7  sin ‘1>) 

+  sin  03S'  sin  03S'  (-  cos  <l>)  +  cos  0gs'  (-  sin  7  sin  <I>), 

and  where  subscript  S'  refers  to  the  illuminated  point  on  the  mean  surface. 

For  some  purposes  (e.g.,  characterizing  target  position),  it  is  desirable  to  have  direct 
transformations  between  the  S'  and  Sz'  systems,  which  do  not  involve  the  angles  5  or  ‘l>. 

To  develop  these  transformations,  we  first  invoke  (1. 10)  and  (I.15.a),  from  which  we  obtain 

*3'  =  (*’  ~  xr  )  cos  y  +  (z1  -  zR ')  sin  7, 

yz'  =  ~  (/-  yR'),and  (I.29.a) 

a3'  =  (x'  -  xR')  sin  7  -  (2'  -  zR ')  cos  7. 
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From  (1. 11)  and  (I.15.b), 

x'  =  xR'  +  jc3'  cos  7  +  z3'  sin  7, 

y'  =  yR'~  y3',  and  (I.29.b) 

z  =  zR  +  jc3'  sin  7  -  zg'  cos  7. 

The  transformations  of  unit  basis  vectors  are  obtained  from  (1.14)  and  (I.16.a,b)  and 
are  as  follows: 


ix'=  ix  '  cos  7  +  X2  'sin  7, 


*  *  —  _  »  f 
iy  “  iy3 ’ 


4*' =  ix3'siny- Xz3' cos^ 


(I.30.a) 


and 

lx  '  =  4x'cos7+  iz'  sin  7, 

3 

iy3'  =  -4y',and  (I.30.b) 

i  '  =  4/ sin  7-  ^'cos7. 

w 

Transformations  between  the  spherical  coordinates  of  the  S3'  system  and  the 
rectangular  coordinates  of  the  S'  system  are  obtained  from  (I.29.a,b)  through  the  standard 
rectangular-to-spherical  relationships.  From  (I.29.b) 

x  =  xR  +  r3'  (sin  0g'  cos  <t>3  cos  7  +  cos  0g'sin  7), 

y'  =  yR-  r3'  sin  0g'  sin  03',  (1.31) 

z'  =  zfl'  +  r3'  (sin  0g'  cos  0g'  sin  7  -  cos  0g'  cos  7), 

and  from  (I.29.a) 

rg'  sin  0g'  cos  03'  =  (*'  -  xr')  cos  7  +  (z'  ~  zR')  sin  7,  (I.32.a) 

r3'  sin  0 g '  sin  0g'  =-  (y'-  yfl'),  and  (I.32.b) 

r3'  cos  03'  =  (jc/  -  xR)  sin  7-  (z  -  zR)  cos  7.  (I.32.c) 
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Substituting  r3'  obtained  from  (I.32.c)  into  (I.32.a,b)  we  obtain 

tan  03'  cos  <t>3'[(x'  -  xR ')  sin  y  -  (z'  -  zR')  cos  7]  = 
(x1  -  xR')  cos  7  +  (2'  -  zR ')  sin  7,  and 


(I.33.a) 


tan  03'  sin  <t>3'[(x'  -  xR')  sin  7  -  (2'  -  zR ')  cos  7]  =  -  ( y '  -  yR')-  (I.33.b) 

humming  and  squaring  (I.33.a)  and  (I.33.b),  then  taking  the  square  root  and  the 
gent  of  the  result,  we  have 
/ 


<t>3'  =  tan'1 


1  of  (1.32. 


.c)  for  r3'  gives  us 


— - ^ -  •  (1.34  .a) 

:R')  cos  7]  |  2 

(I.33.a),  then  taking  the  arctangent,  we  obtain 

_ ~  (y'-yj?') _ 

(x'  -  xR ')  cos  7  +  (z  -  zR ')  sin  7 


xR)  sin  7-  (2'-  2 

;  ratio  of  (I.33.b)  to 


(I.34.b) 


(*'-  x/)sin  7-  (z Zr) 

1  COS0g' 


cos  7 


(I.34.C) 

~  cos  a3 

to  have  the  x'  and  y'  coordinates  of  a  point  S  on  the  surface  (denoted  by 
eras  of  the  spherical  angles  03'  and  <p3'  of  that  point  (denoted  by  03S', 
dar  altitude  (zR '  -  zs')^hR. 


s  relationship  from  (1.31)  and  (I.34.c)  as  follows: 

hR  cos  7  [tan  03S'  cos  <p3S' 

=  jcp'+ - - - 

H  [1  -  sin  7  (tan  0? 


Ks'  cos  7  f  sin  7] 

_ ,  and  (1.35.: 

'3S'  cos  03S'  cos  7  +  sin  7)] 


•a) 


hR  cos  7  tan  03s'  sin  </>3S' 

ys  [1  -  sin  7  (tan  Ooc' cos  cos  7  +  sin  7)] 


(I.35.b) 
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Appendix  II 

SMALL-SCALE  SURFACE  FLUCTUATIONS 

The  averages  over  small-scale  surface  fluctuations  are  calculated  in  this  appendix.  We 
note  in  the  second  section  that  the  ripple  surface  is  given  by  z(g  ,t),  where  =  (x,y)  = 
position  coordinates  on  the  mean  surface,  considered  locally  as  planar 

The  triple  Fourier  transform  of  z(p,f )  is 


Z(n)=ff  (II.l) 


where  12  =  (fe.oj)  =  ( kxky,u ). 

We  assume  that  z(j}  ,t)  is  a  sample  function  of  a  zero-mean,  statistically  homogeneous*, 
wide-sense  stationary  Gaussian  random  process.  From  the  “zero-mean”  assumption  it 
follows  that  the  ensemble  average  (denoted  by  <  >)  of  z(£,t)  is 

<z(g,f)>  =  0.  (II.2) 

From  the  assumptions  of  statistical  homogeneity  and  stationarity,  it  follows  that 

<z(p  ,t)z(p  +  A p,t  +  At)>  =  R(Ap ,At),  both  independent  of  p  and  t.  (II.3) 
Using  (II.2)  and  (II.3)  in  (II.l),  we  have 


oo  oo 


<Z(12)>=//  d^J'  dt<z{^,t)>e-Jut~  mM  =  0,  and  (II.4) 


—  OO  -  oo 


♦An  alternative  terminology  might  be  “spatial’ stationary.” 
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oo  oo  oo 


/*>  JJdSi  2  J'  dt2<^z^l't0z^2't2^> 


-OO  -OO  -OO  -OO 


-'«Vi  +a>2f2  ♦  ft,  *5l  +fe2  *g2)_ 


=  yy*  dg1dt1d(A<g)rf(At)fi(Ag.At) 


(II.5) 


e-/[(u,  ♦  w2)tj  +  (*,+k2)  ‘  Px]  e-;[w2A<  +  k2  •  Ap] 


We  note  that  (See  (2.7.a,b)) 


J f  +  =(27r)26(fclx+fe2x)6(fcly +fe2y),and  (II.6.a) 


yy*  d(j  e  J<W‘  »  2ir  Stciij  +  w2). 


(II.6.b) 


Using  (II.6.a,b)  in  (II.5),  we  have 


<Z(«l)Z(n2)>' 


OO  OO 


{2ir)*ff  di'f  ^(jgvy<Wi* * *£  * 


—  oo  — oo 


where 


6(S2l  +fi2)  =  5(felx  +fe2x)  5<fely  +  ft2y)  5(wl  +  W2^ 


(H.7) 


The  integral  in  (II.7)  is  the  power-spectrum  (in  both  spatial  wave  number  k  and 
frequency  w)  of  z(]),t),  which  may  be  denoted  by  W(S2)  or  W(k  ,co);  thus, 

z(Qi)Z(Q2)  =  w(j2i)5(ni  +n2), 

where  Jli  =  (fti.oji)  and  W(Sl)  has  the  property  W(fl)  =  W(-O). 

k  k  “T  — ^ 


(II.8) 
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Because  z(p  ,t)  must  be  real,  it  is  evident  from  (li.l)  that 


Z*(fi) 


oo  oo 

■JJ-iJ 


dtz(p,t)e 


jut  +  jk  ■  p  _  7, 


Z(-£2). 


(11-9) 


As  a  consequence  of  (II. 8)  ind  (II.9), 

<z(n1)z*(n2,>  =  <z(Q1)Z(-n2)>  =  5^  -  g2).; 

A  general  property  of  zero-mean  Gaussian  random  functions  Xx  2  3  4  is 

<xxx2x3x4>  =  <xxx2>  <x3xA> 

+  <XxX3  >  <X2X4  >  +  <XxXA  >  <x2xz>. 


(11.10) 


(11.11) 


From  (II. 11)  and  the  fact  that  Z(Q)  is  a  zero-mean  Gaussian  random  function  since 
z(p,t),  its  Fourier  transform,  is  such  a  function  and  it  is  a  we'l  known  property  of  Gaussian 
random  functions  that  their  Fourier  transforms  are  also  Gaussian,  we  obtain 


<7(P1)Z(P2)Z(P3)Z(P4)>  =  <z(P1)z(n2)><z(P3)Z(P4)> 

+  <'Z(P1)Z(§3)>  <Z(P2)Z(P4)>  +  <Z(PX)Z(P4)>  <Z(P2)Z(P3) 

Another  property  of  a  zero-mean  Gaussian  random  function  is  that  its  odd-order 
averages  vanish;  consequently, 


(11.12) 


<Z(P1)Z(P2)Z(P3)>  =  0.  (11.13) 

We  make  the  assumption  that  the  small  scale  fluctuations  on  any  part  of  the  mean 
surface  are  only  very  weakly  correlated  with  small  scale  fluctuations  o”  another  part  of  the 
mean  surface.  We  use  the  extreme  limit  of  this  condition  as  an  approximation  on  integrating 
over  the  mean  surface,  i.e.,  we  assume  that 

<fmaHa)smb^b)>*  0  if  |rfl-_r6|<e 

(11.14) 

=  Oif\Ha-rb\>e, 

where  e  is  a  small  positive  real  number  that  limits  the  correlation  region  to  a  very  small 
patch  such  that  the  tilt  of  the  mean  surface  remains  uniform  throughout  the  patch.  If  e  is 
large  enough  to  allow  correlation  between  two  points  on  parts  of  the  mean  surface  with 
different  tilt  angles,  then  a  condition  that  will  later  follow  from  the  ensemble  averaging 
process  (namely  a  correlation  only  between  two  points  with  the  same  c^g)  could  possibly 
apply  to  two  points  on  widely  separated  parts  of  the  mean  surface  but  with  the  same  C£0 . 
This  would  not  make  sense  physically.  Since  the  scale  of  mean  surface  tilt  is  large  compared 
with  the  scale  of  the  small  scale  fluctuations,  this  assumption  is  easily  justified. 
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Another  assumption  we  make  is  that  the  fluctuation  height  z(  p  ,t)  varies  so  slowly  that 
we  can  neglect  its  time  variations  during  illumination  of  a  patch  by  the  radar,  i.e., 

z(g,t)  =  z(  jo).  (11.15) 

It  follows  from  (II. 1),  (11.15)  and  (2.7.b)  that 

Z(Q)  =  Z(k,u)  =  2n  Z(k)&(u)  (11.16) 


where 


oo 

Z(fe,=  ff  dlzWe~jk'£ 


From  (II.4),  (II.8),  (11.10),  (11.12),  (11.13)  and  (11.16),  we  can  write 


<Z(k)>  =  0,  (11.17) 

<Z(k1)Z(k2)>  =  W(kx)  6 (fej  +  fe2),  (11.18) 

<Z(k1)Z*(k2)>  =  W(k1)S(k1  -  k2),  (11.19) 

<Z(fe1)Z(fe2)Z(fe3)>  =  0,and  (11.20) 


<Z( k j )Z( * 2 )Z( k 3 )Z( k 4 ) >  =  <Z(  ft  )Z(  ft 2 ) ><Z( ft 3 )Z( ft 4 ) > 

-2  3  -4  -1  -2  3  -4  (II.21) 

+  <Z(ft  j  )Z(fc  3 ) ><Z( ft  2 )Z(fe  4  )>  +  <Z( ft  j  )Z(* 4 )><Z(*  2 )Z(fe  3  )> 

where  IF(  k )  is  the  spatial  spectrum  of  the  small  fluctuations  on  the  surface  and  each  of  the 
two-fold  product  averages  in  (11.21)  may  be  obtained  from  (H.18|  or  from  (11.19)  and 
(II.9)  when  appropriate.  More  generally,  if  complex  conjugates  Z*(k)  are  involved  in  any 
of  these  averages,  then  (II.9)  can  be  invoked  to  convert  Z*(  k )  into  Z(-k)  and  then  any  of 
the  rules  given  by  (11.17)  through  (11.21)  can  be  applied. 

The  ensemble  averages  needed  in  the  text  can  be  obtained  by  straight-forward 
manipulation  from  (3.7.0,. ..4)  and  (II. 9,  14, 17, 18, 19,  20,  21).  Note  that  those  of  (4.2) 
or  (4.5)  are  all  either  averages  of  a  voltage,  such  as  <UAa2>  and  <UAa4>  (or  complex 
conjugates  of  these)  or  two-fold  product  averages  involving  a  voltage  and  its  complex 
conjugate,  namely  <UAaiU*Bb2'>  anc*  Aal  ^*B63->-  The  first  two  of  these,  by 
virtue  of  (II.8)  and  (11.21),  contain  a  factor  5(a0)  within  their  integrands  that  limits 
contributi-  r  to  points  of  normal  incidence,  i.e.,  those  points  for  which  a0  =  0;  thus,  they 
become  summations  of  contributions  from  such  points  as  given  in  (4.5.3)'  and  (4.5.4)'.  The 
two-fold  product  averages  <U Aa%U* Bi}m>  by  virtue  of  (3.7.1,  2,  3),  each  contain  a 


63 


RAEMER 


factor  FAa(0  3'  ,<l>3',t)FBb(0  3'  >(j)3,t),  which  in  turn  by  virtue  of  (3.5)  and  (11.16)  contains  a 
phase  factor 


±2y*-[V03’*3*'>-V°3  -*3 


I  *  At)) 


Another  type  of  two-fold  product  average,  not  appearing  in  (4.2)  or  (4.5)  but  appearing 
only  in  (4.4.a,...,e)'  and  conspicuou -,ly  absent  from  (4.6),  is  <UAavUBbm>  or  its  complex 
conjugate;  specific  examples  being  as  in  (4.4.a)',  <UAclUBbi>  as  in 

(4.4.b)',  <UAa3UBbl>  as  in  (4.4.c)'  and  <UAa2UBb2>  as  in  (4.4.d)\ This  tjpe  of, 
average,  again  by  virtue  of  (3.5)  and  (11.16),  contains  a  factor  FAa(03',<p3',t)FB b(03 ',<t>3\ 
t  +  At)  implying  a  phase  factor 


±2y  — 
e  c 


*  r  (0  t  +  At)] 


Finally,  there  are  the  three-fold  product  averages  of  the  type 

<VAaAblUCc2> 

2  1 

and  their  complex  conjugates;  specific  examples  being  <U%alUBblUBb2>  and 

as  given  in  (4.4.b)',  <U%blU AolU Aa2>  and  <U*Aal  U%bl  UAa2> 
as  in  (4.4.c)',  and  <UBblU%blU%a2>  and  <UAalU*Aal  U%b2>  as  in  (4.4.d)'.  These  , 
averages,  again  because  of  (3.5)  and  (11.16)  contain  a  factor  FAa(03',tl>3',tl)FBb(Q3  ,<P3  ,t2 ) 
FCc(0 3  ',<j)3',t3),  where  two  of  tiie  tk’ s  inside  the  arguments  of  these  functions  are  equal 
while  the  other  is  in  general  different;  e.g.,  tx  =  t,  t2  =  t,  t3  =  t  +  At.  The  phase  factor,  then, 
has  the  form 


t2j 


lr0^3  ’^3  r0^3  ’^3  ,l2^  +  r0^3  ’^3  ,f3  ^ 


The  eventual  neglect  (in  (4.6))  of  the  averages  <UAa%UBbm>  and  all  terms  of  the 
three-fold  product  averages  except  the  term  of  the  form  <UAaiU%bi><UCc2>  is  based  on 
the  nature  of  the  phase  factors  mentioned  above.  The  time  separation  At  is  sufficiently 
small  to  allow  the  approximation 


±2y~  lro(0g\0g\O  -  r q(8 3' ,<t> ,t  +  A‘0] 


(11.22) 


Note  that  (except  when  At  =  0)  (11.22)  is  not  meant  to  imply  that  the  phase  factor  is 
actually  or  even  approximately  replaceable  by  unity;  it  implies  only  that  it  is  not  so 
oscillatory  as  to  “wipe  out”  the  {d3’,<f>3')  integral.  On  the  other  hand,  the  phase  factors 
that  appear  in  the  averages  <UAaiUBb\>  and  the  three-fold  product  averages  are  much 
more  oscillatory  because  the  values  of  r0  in  the  exponents  (which  of  course  must  be 
positive)  are  added  rather  than  subtracted;  hence,  the  integrands  in  these  averages  are 
products  of  relatively  smooth  functions  and  highly  oscillatory  phase  factors  and  should 
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Ik*  extremely  small.  Based  on  these  qualitative  arguments,  the  following  approximations 
will  he  made: 


<UAaHUl.hrn>*0fat^nm1^ 

<lJAa\U KbitJCc2>  ~  <'^AalUHb}><'t~’cc2'>' 
(with  the  aid  of  (11.21)) 

<'UAa\UliblUCc2>  %  °» 


(II 

(1 


(II 


(with  the  aid  of  (11.21)  and  (11.23)). 


A  final  average  that  appears  in  (4.4.d)  is  (from  11.12)  with  the  aid  of  (11.23)), 


Aa\V AalU Bbl^BhV*  <UAal  UAal><UBb  1  U!ib  1 


(II 


Appendix  III 
THE  MEAN  SURFACE 

The  mean  surface  or  “swell”  surface  is  characterized  by  specifying  angles  6  and  <I>  as 
functions  of  03'  and  03\  We  will  first  specify  the  mean  surface  as  a  function  of  x’  and  y', 
denoted  by  H'(x',y'). 

Referring  to  Appendix  I,  we  note  that  a  point  on  the  mean  surface  must  by  definition 
fulfill  the  condition  2  =  0,  because  the  mean  surface  is  explicitly  defined  as  the  z-  0 
surface.  We  set  2  =  0  in  (1.6)  and  note  that  2'  in  (1.6)  is  functionally  related  to  (x',y') 
through  the  equation 


z  =  H'(x',y')  for  2  =  0. 

Referring  to  Figure  (1-2)  and  Eqs,  (1.1)  and  (1.2),  we  note  that 

x”  =  (x'  -  Xq  ')  cos  <l>  +  (/  -  yQ')  sin  <I> 
/"  =  -  (x'  -  x0')  sin  d»  +  (/  -  yQ')  cos  4> 


(III.l) 


(IH.2) 


Representing  x'"  and  y'"  as  functions  of  x'  and  y'  through  (III.2),  writing  H'(x',y')  in 
the  form  H'(x"'(x',y'),y"'(x,,y')),  and  applying  the  chain  rule,  we  obtain 


But  from  (III.2) 


dH'  dH'  dx'"  dH'  by'” 

-  —  -  - - -  4.  -  ■■■■■  — —  onr] 

dx'  dx'”  dx'  dy'"  dx'  ’ 

dH'  dH'  dx'”  dH'  dy'" 

—  -  + 


dy'  dx'"  dy'  dy'”  dy'  ' 


•\  fff 

dx  a 

— -  =  COS  <P 

dx 
dx'" 


dx' 


=  sin  4> 


9y‘ 


dx 


—  =  -  sin  (1>,  and 


3y 


9y 


—  =  cos  d>. 


(III.3.a) 

(III.3.b) 


(III.4.a) 

(III.4.b) 

(III.4.C) 

(III.4.d) 
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Since  x"  is  the  direction  of  the  local  tilt  of  the  surface  relative  to  the  horizontal 
and  y'"  is  the  direction  normal  to  that  tilt,  it  follows  that 

bH' 

~~177  ~  tan  6 ,  and  (lll.S.a) 

dx 


a//' 


bv' 


=  0. 


It  follows  from  (III.3.a,b),  (III.4.a,...,d)  and  (III.5.a,b)  that 

bH 

H  ,  =  -  tan  5  cos  <1>.  and 

x  bx 

bH 

H  ,  =  —  =  tan  5  sin  <1>. 

•v  ay 

Squaring  (III.6.a,b),  adding  and  taking  the  square  root,  we  have 

tan  5  =  |v'H  |  =  y/‘(Hx,)2  +  (Hyl) 2, 

or  equivalently, 

1 


cos  §  =  - 


sifi  5  = 


VT+  IV'H  K 
IV'// 1 


,  and 


y/i  +  lv'//l2 


(U1.5.b) 


(Ill.e.a) 

(in.e.b) 

nn.7) 

(III .8  .a) 

(III.8.b) 


Dividing  (III.6.b)  by  (III.6.a), 


or  equivalently, 


tan  <1>  = 


.  "«• 

cos""iv'w'and 


sin  <l>  = 


IV'// 1  ' 


(III.9) 


(Iil.lO.a) 


(Ill.lO.b) 
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We  choose  as  a  model  for  H(x',y')  a  sinusoidal  function  of  the  form 

H(x',y')  =  C  cos  (Kxx'  +  K2y'  +  *).  (111.11) 

Differentiating  (III.ll)  with  respect  to  x'  and  y',  we  have 

Hx,  =  -  KXC  sin  (Kxx  +  K2y'  +  <//),  and  (111.12*) 

Hy '  =  -  K2C  sin  (Kxx'  +  K2y'  +  i//).  (III.12.b) 

Using  (111.12.a,b)  in  (III.8.a,b)  and  (Ill.lO.a.b), 

1 

COS  8  =  == . . .  ,  (IH.13.a) 

V  1  +  C2(K2  +  K\)  Sin2(ff1*'  +  K2y'  +  </,) 


C  s/K\  +  K\  sin  ( K.x '  +  K2y'  +  4>) 

sin  6  =  -==r- ,  (III.13.b) 

V  1  +  C2(K2  +  K\)  sin2(/f1x'  +  K2y'  +  tfr) 

cos  4>  =  ,  and  (III.14.a) 

VK2  +/f2 

-JC2 

sin  <1>  =  —=======  .  (III.14.b) 

%/«!♦«! 

To  obtain  cos  5,  sin  5,  cos  ll>  and  sin  <t>  as  functions  of  03'  and  $3'  for  fixed  values 
of  7  and  h^,  we  use  (I.35.a,b)  in  (III.6.a,b)  and  (IH.lO.a.b).  If  we  employ  the  specialized 
model  (III.ll),  then  cos  4>  and  sin  as  given  by  (III.14.a,b)  are  independent  of  x'  and  y' 
and  cos  8  and  sin  5  can  be  expressed  as  functions  of  03'  and  03'  by  substituting  (I.35.a,b) 
into  (III.13.a,b).  Summarizing  these  results,  first  for  the  general  case,  we  write:  (from 
(III.8.a,b)  and  (III.10.a,b)) 


cos  5  =  — 7=  =•  , 

Vi  +  Iv'hI2 

„  \v'h  i 

sin  8  =  — .  —  , 

Vi  + Iv'hI2 


(III.15.a) 

(III.15.b) 
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COS  <1*  = 


IV'// I 


,  and 


sin  <l>  = 


//  , 

y 


IV'// 1 

where  (from  (I.35.a,b),  with  removal  of  the  subscripts) 


(II1.15.C) 


(II1.15.d) 


3// 


a// 


//  ,  =  ,  H.  =  , 

-v  3y 

hR  cos  7  (tan  03'  cos  03'  cos  7  +  sin  7) 

x'  =  x  o'  + - ; - ; - 

n  [1  -  sin  7  (tan  03  cos  03  cos  7  + sin  7)] 

hR  cos  7  tan  03'  sin  03' 


(111.15.e) 


(III. 15. f) 


y  =yR 


[1  -  sin  7  (tan  03'  cos  03'  cos  7  +  sin  7)] 


and  H(x',y')  is  an  assigned  real  function  of  x'  and  y'.  For  the  case  where  (111.11)  is  used 
(from  (III.13.a,b)  and  (III.14.a,b)), 


cos  5  ~ 


sin  5  = 


y/l  +C2(K*  +  K\)  sin2  (ATjjc'  +  K^y'  +  0)  ’ 

Cy/K\+K\  sin  ( Kxx '  +  K2y'  +  0) 

sin2  (K jjc'  +  K2 y'  +  0)  ’ 

-  /c. 


cos  <i>  = 


vA? 


,  and 


+  Ki 


~Kn 


sin  <l>  = 


V k\+k\ 


(III.16.a) 


(III.16.b) 


(1II.16.C) 


(I11.16.d) 


where  x'  and  y'  are  given  in  (III.15.f). 
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Some  of  the  terms  in  our  final  results,  (UAa 0,  <UAa2>  and  <UAa4>)  contain  the 
condition  a0  =  (aox  ,aoy )  =  0.  From  (1.28)  (dropping  subscripts)  this  condition  is 
equivalent  to 


tan  0 3  cos  03'  (cos  5  cos  7  cos  «t>  -  sin  5  sin  7)  +  tan  03'  sin  03'  (-  cos  6  sin  <1>) 


=  -  (cos  5  sin  7  cos  4>  +  sin  5  cos  7), 


(III.17.a) 


tan  03'  cos  03'  (cos  7  sin  $)  +  tan  03'  sin  <p3  (cos  4>)  =  -  sin  7  sin  <J>,  (II1.17.b.) 


tan  63'  cos  03' 


(sin  5  cos  7  cos  <fr  +  cos  6  sin  7) 
(cos  5  cos  7  -  sin  6  sin  7  cos  4>) 


(tan  3  cos  $  +  tan  7) 
(1  -  tan  5  tan  7  cos  «l>) 


,  and 


tan  03'  sin  <f>3 


_ (sin  6  sin  4>) _ 

(cos  6  cos  7  -  sin  5  sin  7  cos  4>) 


(tar.  6  sin  $  sec  7) 

(1  -  tan  5  tan  7  cos  «l>) 


(III.18.a) 


(III.18.b) 


Squaring  and  summing  (III.18.a,b)  and  then  taking  the  ratio  of  (III.18.b.)  to  (III.18.a) 
we  obtain 


(0a\  =  o  =  ten 


-1  f  V  (k*0  &  cos  y  cos  ^  +  s*n  7)2  +  (fan  B  sin  4>)2 


(Odo  =  0  =  tan'1  | 


(cos  7  -  tan  5  sin  7  cos  4>) 
-  tan  5  sin 


,  and  (111.19*) 


(III.19.b) 


tan  6  cos  7  cos  «I>  +  sin  7 
If  the  surface  is  perfectly  horizontal,  i.e.,  5=0,  then  Eqs.  (III.19.a,b)  degenerate  into 


(03')d  s  0  =  tan"1  I  v/^-7)2  |  =  7,  and  (III.20.a) 

0  u  (  cos  7  ) 

5  =  0 

(<t>3%  =  0  =  tan"1  j -  ~  }  = (III.20.b) 
5  =  0  ' 

The  implication  of  (III.20.a,b)  is  that  the  angle  03'  corresponds  to  the  vertical  direc¬ 
tion  (see  (I.15.a,b)  and  Fig.  1-5).  Thus  makes  physical  sense  because  the  condition  %0  =  0 
corresponds  to  normal  incidence  on  the  surface,  which  would  obviously  occur  only  with 
vertical  incidence  if  the  surface  were  perfectly  horizontal. 
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The  condition  on  the  angles  (03', tp3')  corresponding  to  a()  =  0  can  be  determined 
from  (III.17.a,b)  or  (III.18.a,b)  or  (III.19.a,b)  combined  with  (III.15.a,...,f).  The  results  in 
the  general  case  are:  from  (III.18.a,b)  and  (III.15.a,...,d) 


-  (Hx,  +tan7)a-(j  „  0 

(ten  d~'  cos  <bJ)  *  -n  =  - — - ,  and 

3  '«o-°  (1-//  ,tan7)^ 


a  =0 


(II  I.21.a) 


(Hy')a  =  osec? 


(tan03'sin03')^  .  0  « 


(1  -  Hx,  tan  7)5  =0 


III. 21. b) 


where  (from  (III.15.f))  Hx<  and  Hy  >  are  functions  of  the  two  variables 

h R  cos  7  (tan  0„ '  cos  03 '  cos  7  +  sin  7)^  _  0 

(x ')'  .  0  =x'R  + - - - —7 - -1 - ,  and  (III.21.C) 

So  0  R  [1  -  sin  7  (tan  03  cos03  cos  7  +  sin  7)]  -  _0 

-*o 

cos  7  (tan  03' sin  ^>3')^  =0 

(y')d  ,0=y'R-  7;— — 77— rr - tt - - •  (in.2i.d) 

~o  w  n  [1  -  sin  7  (tan  03  cos  03  cos  7  +  sin  7)]  ^  _  0 


In  the  special  case  where  the  model  (III. 11)  is  used,  with  the  aid  of  (III.16.a,...,d), 
Eqs.  (IIT.21.a,b)  have  the  form 


{CKX  sin  (K xjc'  +  K2y'  +  0)  -  tan  y)&  =  0 


^an  °3  cos  °3  )af>  =  0  [1  +  CKX  sin  (Kxx'  +  K2y'  +  0)  tan  7]  ^  =0 


,  and  (III.22.a) 


-  CK2  [sin  (Kxx  +  K2y' +  =  0  sec  7 

t  I  ~*n 

(tan  0 o  sin  0,  )  -  _  n  = - ; - ; - , 

3  ^3  -  0  [i  +  CKX  sin  (/C^  +  K2y  +  0)  tan  7]  -  =0 


(III.22.b) 


where  x'  and  y'  are  given  by  (III.21.c,d), 

In  describing  the  position  of  a  point  on  the  surface  illuminated  by  the  radar,  one  can 
consider  this  surface  as  perfectly  horizontal  except  where  phase  differences  are  involved. 
The  quantity  (zR  ’  -  z')  has  usually  been  considered  equivalent  to  hR,  the  radar  altitude, 
because  the  difference  is  usually  negligible.  If  the  highest  swell  were  30  m  (100  ft)  and  the 
radar  altitude  were  1.6  km  (1  mile),  the  error  in  considering  ( zR '  -  z)  equivalent  to  hR 
would  be  about  2  percent.  In  most  cases  the  error  would  be  even  smaller. 

The  phase  factor  e/fer«  which  appears  in  the  field  is  “washed  out”  in  all  of  our 
averages  exept  UAao,  <UAa2>,  and  <UAa 4>.  Thus,  we  need  only  consider  the  effect  of 
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this  factor  to  calculate  those  terms.  Referring  to  the  diagram  of  Fig.  III-l,  the  vector  rt0  is 
that  from  the  illuminated  point  [*',  y',  on  the  mean  surface  to  the  radar.  The 

vector  40  is  that  from  the  point  (x',y\0)  (directly  below  the  illuminated  point)  to  the 
radar.  The  vector  f}'  is  equal  to  _(.2'//'(x',y’),  where  \z'  is  the  unit  basis  vector  in  the  z 
direction.  It  is  evident  that 

from  which  it  follows  that 


where  0R'  is  the  polar  angle  of  the  radar  with  respect  to  the  point  (x',y',0),  described  in 
the  {x',y'jz')  coordinate  system. 


Fig.  III-l  —  Mean  surface  geometry 
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To  express  the  phase  factor  entirely  in  terms  of  the  coordinates  we  note 

that  (with  the  aid  of  (I.35.a,b)  removing  the  subscript  S ) 


(111.25) 


r0  =y(xR'-x')2HyR'-y')2+hl 

hR  cos  7  sec  0 3' 

(1  -  sin  7  (tan  03' cos  03' cos  7  +  sin  7)] 

It  follows  from  (III. 24)  and  (III. 25)  that 

hR  cos  7  sec  03' 

r  - - - - - 

°  [1  -  sin  7  (tan  03'  cos #3'  cos  7  +  sin  7)]  (III  26) 

-  sec  7  cos  03'  [1  -  sin  7  (tan  03'  cos  03'  cos  7  +  sin  y)]H(x\y'), 

where  x  and  y'  in  H(x',y')  are  given  by  (III.15.f). 

The  factor  e*k r<>  jr0  which  appears  in  our  averages  can  be  approximated  through  'III. 26) 


by 


ikr„ 


( 


jkhR  cos  7  sec  03' 


[  1  -  sin  7  (tan  d3'  cos  <j>3  cos  7 


+  sin  7)1  ^ 


rQ  hR  cos  7  sec  03' 

[1  -  sin  7  (tan  03'  cos  03'  cos  7  +  sin  7)] 

-jk  sec  7  cos  0  J  [1  -  sin  7  (tan  0 cos  <j>J  cos  7  +  sin  y)]H(x',y') 

fi  «  OO 


(III.27) 


where 


*  =XR  +  ■ 


hR  cos  7  (tan  03'  cos  03'  cos  7  +  sin  7) 


[l  -  sin  7  (tan  03 '  cos  <j>3 '  cos  7  +  sin  7)] 
hR  cos  7  tan  03'  sin  03' 

yf  =  y  ^  „  -  -  - ■■  —  ■  —  1  ■-  1  -  — 

R  [1  -  sin  7  (tan  0g' cos  </>3' cos  7  +  sin  7)] 
(again  with  the  aid  of  (I.35.a,b)). 


,  and 
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In  the  special  case  where  (III. 11)  is  used,  (III.27)  has  the  form 

/  jkhR  cos  7  sec  03' 


jkr  \  f  1  -  sin  7  (tan  7  ’  cos  0,’ cos  7  +  sin  7)] 


rQ  hR  cos  7  sec  03' 

[1  -  sin  7  (tan  03'  cos  03'  cos  7  +  sin  7)] . 

-jk  sec  7  cos  O'  (1  -  sin  7  (tan  0  *  cos  0,'  cos  7  +  sin  7)] 

9  d  3d 


(III.28) 


Ceos  <  K2yR'  + 

cos  7  (tan  03'  cos03'  cos  7  +  sin  7)  -  K2  tan  03'  sin  03'] 

[1  -  sin  7  (tan  03 '  cos  03 '  cos  7  +  sin  7)] 

In  the  terms  where  the  factor  ejk  r°  /r0  is  not  washed  out  in  the  averaging  process, 
i.e.,  UAa,  <UAa2>  and  <UAai>,  the  condition  a0  =  0  prevails.  Consequently,  when 
(111.2*7)  is  used  in  these  terms,  the  angles  03'  and  03'  obey  the  conditions  (III.21.a,b,c). 
When  (III.28)  is  used  in  these  same  terms,  (i.e.,  when  the  model  (III. 11)  is  applied)  then 
03'  and  03'  obey  the  conditions  (III.22.a,b)  and  (III.21.c,d). 


Appendix  IV 

CALCULATION  OF  FIRST-ORDER  RESULTS 

To  obtain  the  result  (5.1),  we  first  specialize  to  the  case  where z(p,t)  does  not  vary 
significantly  with  time  during  the  period  of  illumination;  hence,  we  can  denote  z{j),t)  by 
z(p)  and  its  Fourier  transform  by  Z(i$),  having  eliminated  by  this  specialization  the 
necessity  of  Fourier  transforming  with  respect  to  time  as  well  as  the  spatial  variables. 

From  (2.16.1)  (with  the  assumption  indicated  above) 


#<1>(4,0J)  =  J j' dgd03dj3l  dtJC 


jc  *  j?l  "  e(Wl  "  1 


(IV.l) 


*(  1  i  ^ )  !/(|t  ^ 


where 


fi(1)(|,w/j§1,w1^1,f1)- - 7 - j  T  &(£)a(j3,fe  ), 

(2ir)3  A(  Ig  I) 


and  where  b(p)  is  the  4  X  4  matrix  given  in  (2.16.0)';  A(  l£  I)  =  (j32  +  Yz)  (1  -  /jf  +  (lz  Yz), 
and  a(P,k  x )  is  a  4  X  2  matrix  defined  by 


6(0,) 


(IV.2) 


where  B0(^x)  is  obtained  from  (2.12.a)"  and  (2.12.b)"  with  J3  set  toJJj,  and  the  elements 
of  the  4  X  4  matrix  A1(j5,k1)  and  the  4  X  2  matrix  &1(Q,k1)  are  obtained  from  (2.12.a)"  ’ 
and  (2.12.b)"  ’  respectively.  The  latter  elements  are  the  first  order  terms  of  the  series 
expansion  of  (2.12.a)" '  and  (2.12.b)" '.  The  first  order  terms  of  (2.12.a)'' '  are  as  follows: 


(IV.3.a) 


<412(1))l  =J2  Zy->?PzZ 

-7* 

'z 


(IV.3.b) 

(IV.3.C) 
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<V’>. 

r,  >  w 

(IV.3.d) 

Ms.'1’*! 

-  Py2y  ~  if  *<*  -  O 

(IV.3.e) 

(A32^1’)i 

=  “  ^[Zy  +j£PyZ] 

(IV.3.f) 

^33* 1 

=  -  0yzy  +j?z(v2  -  Py2) 

(IV.3.g) 

=  PXlZy  +jfzpy) 

(IV.3.h) 

The  first  order  terms  of  (2.12.b)" '  are: 

<»n(1|)i 

K 

mv> 

(IV.4.a) 

(Biz"1), 

h 

~yzy~  jf^zZ 

(IV.4.b) 

•  ~PyZy  +irza-py2) 

(IV.4.C) 

(B32(1|)l 

=  PjZy  +j?PyZ] 

(IV.4.d) 

To  obtain  the  first  order  results,  we  follow  the  analysis  using  (2.15.0),  (2.15.1), 
(2.16.1),  (2.17),  (2.18.a,b)  and  (2.19)  with  the  aid  of  (2.16.0)':.  Expressing  z(g), 
zx(q)  and Zy(j) )  in  terms  of  its  Fourier  transform,  we  have 


z(/0)  =  f f  dke>b'RZ(k ), 

(IV.5.a) 

2 *(£)  =  J  J *kjkx  e>b  '  £  Z(k),  and 

(IV.5.b) 

Zy  (P  )  =  f  f  d&ky  eJ*  zty- 

(IV.5.C) 
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From  (IV.5.a,b,c)  and  (IV.3.a)  the  element  A( 1  >  could  be  expressed  as 


j  f  fdke>k.  'Bjkx  Z(U) 
ffdk  e'*  '  H  Z(k)(Au{^,k))v 


( IV  .6) 


where 

K 

(An(^i,fe))i  =  jky—  . 

The  elements  of  (Ajk) j  and  ( Bjk)l  obtained  from  (IV.3.a,...,h)  and  (IV.4.a,...,d)  as 
indicated  in  (IV.6),  are  as  follows: 


(An(£>hi))i 

"  $xhly&,~X 

(lV.3.a)’ 

(A12<M  1»1 

(IV.3.b)' 

=  jpx*lyy2~1 

(IV.3.C)' 

(^14(j5’*i))i 

mtfykiy  y/'-ify. 

(IV  .3  .d)' 

(^31(|,fcl))l 

=  iPykly-j?a-fiy2) 

(lV.3.e)' 

"  -  >Pxkly  -  jfMy 

(IV.3.f)' 

(A33^£’£l^l 

=  -jPykly+j?(»2-fiy2) 

(IV.3.g)' 

(A34(£’£l^l 

=  ®xkly  +ifWy 

(lV.3.h)' 

(Bu(Ml)h 

=  Vi/;1 

<IV.4.a)' 

(*12(£4l»l 

=  y^-x-ifh 

(IV.4.b)' 

=  -#,*iy+/r<l-0y2> 

(IV.4.c)' 

(B32(Mi»i 

=  Vly  +  ^My 

(IV.4.d)' 
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The  other  elements  of  A  and  B  are  obtained  through  the  relationships  implied  by 
(2.11.a,b)',  i.e., 


A21  ~ 

~  A12’ 

(IV.7.a) 

A22  = 

~A1V 

(IV.7.b) 

A23  ~ 

-  A  ' 

■*14* 

(IV.7.c) 

*** 

A24  ~ 

-A  ' 

*13’ 

(IV.7.d) 

^41  * 

A  ' 

*32* 

(IV.7.e.) 

it 

<N 

IX 

V. 

(IV.7.f) 

A43  “ 

A34* 

(IV.7.g) 

^44  = 

7  ' 

*33’ 

(IV.7Ji) 

II 

H 
„  Cl 

'V  t 

-*12* 

(IV.8.a) 

b22« 

~*1V 

(IV.8.b) 

*41  = 

Bgjj'.and 

(IV.8.c) 

B42  = 

V. 

(IV.8.d) 

where  the  prime  indicates  that  x  and  y  components  are  interchanged  in  the  corresponding 
original  element,  e.g.,  if  Ajk  =  Pxkly,  then  Ajf{  =  0yklx. 

We  now  invoke  (2.17),  (2.18.a,b)  and  (2.19)  with  the  aid  of  which  we  can  use  (IV.l) 
to  express  the  horizontally  and  vertically  polarized  components  of  the  electric  field  of  the 
reflected  wave  to  the  first  order  of  perturbation.  The  result  of  this  operation  is: 

r*rv<S>i. 

H  H 


(IV.9) 


/ 
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where 


Zy  =  transpose  of  , 

H  H 

Q(0)  is  given  in  (2.17),  and  [£(0)]  i  is  the  first  order  term  of  E(j3, oj)  as  given  by  (IV. 1). 
From  (IV.l),  (IV.9)  and  (2.1(f), 

(£rV(|)ll=  / f  d|itSv(|;|1)]1|i(|1)Z[-^(|i-j^)],  (IV.10) 


where 


[sv(p^)}1  =^iTv 

H  H 


(IV.10)' 


and  where  (a)  is  the  quantity  a(0,fe , )  defined  in  (IV.2)., 

^  Examining  the  definition  of  a( J5,k)  in  (IV.2),  we  note  the  matrix  product 
(■4l(Ml)M£l  ))>  which  is^a  4  X  4  matrix.  The  first  step  is^the  calculation  of  that  product. 
Calling  this  matrix  product  d ,  using  c  to  denote  the  matrix  A x ,  and  invoking  (IV.3.a,...,h), 
(IV.4.a,...,d),  (IV.7.a,...,h),  and  (IV.8.a,...,d),  we  obtain 


'11  C12  c13  c14l  f  bll  b12  b 


13  b14 


1  C12  ~C11  C14  -C1 3  ~b'\2  611  b14  bl3  ttdjk]  d 


C31  c32  c33  C34 


bll  b32  b13  b14 


A  ’a  <IVU) 

A  A 


C32  C31  C34  c33  "b32  -bll  b14  b13 


es  c  and  b,  we  have  (Note  that  b11  = 

:  bll»  b14  : 

Id  n 

“  bn<cn 

+  Ci3)-  h12c!2' 

'  b32c14» 

id  12 

bll(c12 +c14^  +  b12cl] 

L  +  b32C13 

'<*13 

=  b13^Cl 1 

+  cl3)  +  bi4(ci2 

+  C14>* 

'<*14 

=  b14^cll 

+  c13)  +  b13(C12 

+  C14>’ 
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d21  ~  d12 > 
d22  ~  dll’ 

d 

d23  =  "  d14’ 
d24  =  "  d13* 


id2i 

=  bll(C31 

+  c33^  ~  ft12C32  “  ^32c34’ 

(IV.12) 

id  32 

-  bll(c32  +  c34)  +  b12czx  +  632^33, 

id33 

613^C31 

f  C33^  +  b14^C32  ^  C34^» 

>d34 

=  ^14^c31 

+  c33^  +  b13^C32  +C34^’ 

d41  =  "  d32’ 
d42  =  "  d3l» 

^43  =  ^34 
d44  =  rf33’ 

Calculating  the  elements  of  the  matrix  d  from  (IV.3.a,...,h),  (IV.1.4.a,...,d),  (IV.7.a,...,h), 
(IV.8.a,...,d)  and  (IV.12),  we  obtain  (explicitly  indicating  the  arguments  of  the  functions): 

d»  '  <!•*> 

*u<M>  -  ‘Wlgi)*/)’  +-7t*ii»<l/tlX?  +  ^uc<l$lH. 

-  ?13.(  IfilX'A  +  7  W  l&lwxA. 

-  ii4.<lsl>‘A  +7  +?i4c<ISI)i. 
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where 


d32</  +  c '^326^  l£l^v  +  ^32(^i£l^  ' 

*3 8<M>  "  <?33.<l£IV,  *  “|rf,l»(ISI»?  *?33,<l£l)l- 

4<M>  -  d34j Ignkjd,  +  jj34t( IU.DMv 
d21  =  d12> 
d22  =  dl 1  * 
d23  =  "  d14  ’ 
d24  =  "  d 1 3’ 
d41  d32' 

d4  2  =  “  d31’ 
d43  =  d34 ’ and 
d44  =  d33’ 


(IV. 13) 


dlla(l£D  =  I"*'2. 

b 

duc(l£l>  -  0. 

<WI£I> 

b 

d130(|£D  "  ^4.(|£l)-fl*,+7.). 

b  b 


<Wl£l>  =  M\£\) 
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=  <W  l&l  - 

b  b 


A(|g|) 

A(|^|)  +  (l-  v2Mz-  jz) 


d32cC0  =  -  v2Wz  +y2), 
y  Jo  ] 
33“~|1-*2  ’ 


d33c  =  0/M1  "  v2)’  80(1 
d’"  “[-2  - 1]‘ 


(IV.13)' 


Denoting  Sx  by  e  and  3a  by  f  and  invoking  (IV.12),  IV.13)  and  (IV.13)',  we  obtain 


^11  ^12 


_e12  ~eM 


e31  e32 


e32  e31 


a=](e  ~  ~~f)~ 
A 


*11 

dl  2 

d13 

d14 

“12 

dil 

~*14 

w' 

d13 

CO 

>-* 

d32 

d33 

d34 

d32 

'V 

_d3 1 

%4 

d33 

(IV.14) 


where 


«ll“J?ll  =  fflz'lkytX’ 

.  Up* 

el2~  ieX2  ~  iPZ  (kyPy  “  )* 


e31=  i? 31  “  ~iPz  lkyPyPz-  —  (1-Py  *)1. 
.  cofl 

e32=  /e32  =  jPz  Px^kyPz+~Pyh 


1  “ 


(IV.14)' 
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Expanding  (IV. 14),  we  have 

Aall  =  Jl  ell^  +  ^12  ~  d  13^31  “  rf14^32^’ 

^  21  =  e  1  2^  +  ^1 1  +  d\  4^31  +  ^1  3^3  2  ^  ’ 

^a31  =  ][e31^  +  ^32~  d33*3\  ”  ^34^32^’ 

^a41  =  Jle32^~  d31  ~  d34^31  "  ^33^32^’ 

(1 

^a12  =  e12^  '  ^1 1  ~  d  13^32  "  ^1  4^31  =  ^a21  ’ 

Aa22  ~  -  di2  +rf14^32  +rf13^31  =~  ^all’ 

^a32  =  e32^  ~  ^31  ~  d33^32  "  d34^ 31  =  Aa41’  ancJ 

Aa42  =  e31^  +d32  ~  d34r32  ~  d33^31  =  Aa31  ■ 

From  (IV. 13),  (IV. 13)'  and  (IV. 15)  we  obtain,  after  considerable  manipulation  (again 
explicitly  indicating  function  arguments). 


(IV.  15) 


"  A(j|j)  V(M>. 


(IV.  -I 


where 


jn<M)  -  sw  mv. +~°ns(  m  »a  • 

a2l(|>*)  =  a21a-l£l^x^x  +'~Ia21&^£l^x2  +  a  21  l£.i^3  ’ 

“31  =  «31a<I^DfeA  +~f®31^l^|)Py2  +  «31c<I£IH* 

a41^’^)=  “41a(|j^l)fex^y  +  7“  41  > 


pga 


a12 

II 

1 

to  " 

M 

t 

u22 

II 

1 

c* 

►-» 

H-* 

a32 

=  a^.and 

a4  2 

=  a31’ 
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and  where 


2\i„<l£l>  =  2(p2-  l), 

b 

«21aHPj)  =  2(^2  -  1), 

b 


6  2(*» 
«aic(lfiD  =  ‘  2(^2-l)72, 


0 


2(c  -  1)(72  -  Pt) 


'  41a 
b 


(!£i)  = 


o 


-  2(^2  -  1)(72  -  Pz) 


a22(M>  = 

®32(f’£)  =  and 

«42<M>  =  («3l(M»'- 

Using  the  above  results  we  can  write  (IV. 16)  more  compactly  as  follows- 

j 


where 


«u<44>  -  ^|3'h.<KK  +  f‘ nudiiw,!- 
-  - 4[ifu(M>i 1 + y  »2ic<iai). 

■^(®4l(M)1'  +T-a3icO£lKand 


-  «,K».<I4I»,  +f ' «3i»<ISlw,  ;■ 


(IV. 16)' 
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To  specialize  the  results  to  the  backscattering  case,  we  first  note  that  the  unit  vector  aQ 
is  defined  as  that  directed  from  the  radar  to  the  illuminated  patch.  In  vector-matrix  notation 


where  oloz  =  -  N/l"~-  aox2  -  aoy2 


a ,  = 

-*o 


a. 


or 


oy 


a 


(IV.18.a) 


We  now  define  a  two-dimensional  unit  vector  a0,  the  projection  of  a0  on  the  x-y 
plane,  as 


Of 

->o 


cl 


oy 


(IV.18.b) 


For  backscattering  of  an  incident  plane  wave  or  plane  wave  pulse,  the  unit  vector 
corresponding  to  the  direction  of  the  incident  wave  is  + 1£0,  while  the  unit  vector 
corresponding  to  the  direction  of  the  backscattered  wave  is  -  o;0  (see  2.1c).  Then, 


h  - 


@lx 

Ply 


a. 


ox. 


Of.. 


<-*y 


=  ao,and 


(IV.19.a) 


/»  = 


L 

J 

L  J 

r  i 

B 

-  O' 

X 

= 

ox 

1 

. 1 

-  rv 

=  -  O' 

-+o 


(IV.19.b) 


Continuing  the  specialization  to  the  backscattering  case,  we  invoke  (IV.19.a,b)  in  the 
expression  (IV. 10)'  resulting  in  an  expression  of  the  form 


2  co  . 


««Hi  "  I»V(-  So>- 


(IV.20) 


H 


H 


To  continue  the  specialization  to  backscattering,  we  invoke  (IV.19.a,b)  in  the 
expression  (IV  .9)  and  carry  out  the  matrix  operations 


s.)U  -  i? «-  »„>M- 


H 


H 


% 

t. 

Jet?'' 


(IV.21) 
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where  B(d0)  is  a  4  X  2  matrix  which  will  be  specified  below,  where  £y  and  E  H  are  given 
by  (I.24.a,b)  or  (1.27),  Q(~  aQ)is  given  in  (2.17)  with 

y  y 

02  =-  aoz,  b0(-  a0)  is  [^(J^))'1  as  given  in  (2.16-0)'  with 

0*  ~  aox> 

y  y 

0 z-~  aoz,  and  L(a0)  is  given  in  (2.19)  with 

0  =  a  ,  0  =  -  a  .• 

HX  OX,KZ  02 

y  y 

•*** 

The  indicated  operations  lead  to  an  expression  of  the  form  (3.7.1).  The  E's  appearing 
in  (IV.9),  (IV.  10),  and  (IV.20)  are  designated  are  U’s  in  (3.7.1)  and  the  quantities 

<v 

Sv,  etc.,  appearing  in  (IV.20)  are  the  quantities  SAal 
H 

as  given  in  (3.7.1),  where  Aa  =  VV,  VH,  HV  or  HH,  depending  on  incident  and  received 
polarizations. 

These  quantities  SA  a  j  are  obtained  by  carrying  out  the  calculations  indicated  in 
(IV.21)  with  the  aid  of  (2.16-0)',  (2.17),  (2.18.a,b)  and  with  the  appropriate  specializa¬ 
tion  of 


0X ,  as  follows: 

y 

2 


tsvv(-“0)l 

VH 

HV 

HH 


± _ 1_ _ 

*  *'<&*&>. 


(Syy) 

~  (Syy, 

VH 

VH 

HV 

HV 

HH 

HH 

(IV.22) 
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where  (primes  on  the  elements  {Bl{)[m  indicate  jc-.v  interchange) 

(•  {a«z  +  WKn(/^11  +a,>x^'>\\ 

1  T  ' 


®vv> 

=  e'' 

kh  X 

,xi'vh 

VH 

V 

V 

HV 

H 

H 

HH 

V 

1 

L 

H 

-  a  , 
oy 

V 

H 
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V' 

// 
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,  and 

r 

V 

V 

H 

H 

> 

H 

H 

>; 

£ 

i 

>• 

i 

( svvY 

~  same  as  (Svv,) 

VH 

VH 

HV 

HV 

sf 


HH  HH 

except  that  x  and  y  are  everywhere  interchanged. 

(a0*J  ')z  =  2  component  of  cross-product  of  a0and  fi/,  equal  to  (a„x£y  ~  %y^x)- 

A'  =  (y02-  ao^(1~aL-  Vo2)  =  A(Hoi) 

(SvX$.hK  =  ^Vx^Hy  ~  ^Vy^Hx* 

The  quantity  [Svv(-  a)]  given  in  (IV.22)  is  the  element  of  the  first  order  “polariza¬ 
tion  matrix”  corresponding  to  vertically  polarized  reception  and  vertically  polarized 
transmission,  where  the  illuminated  surface  patch  is  at  an  angle  a„  with  respect  to  the 
radar  The  subscript  VH  indicates  vertically  polarized  reoeption,norizontally  polarized 
transmission,  HV  indicates  horizontally  polarized  reception,  vertically  polarized  trans¬ 
mission  and  HH  indicates  that  both  reception  and  transmission  are  horizontally  polarized. 

We  now  turn  to  the  evaluation  of  the  matrix  element  B  appearing  in  (1V.22).  In  the 
first  order  case,  these  are  the  elements  of  the  matrix  a  given  in  (IV. 17)  with  k  set  equal  to 
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-  —  a Q  and  ^  set  equal  to  a Q .  From  (IV.  16)',  (IV. 17),  and  (IV.22)  with  the  above 
specifications  of  the  variables,  we  have 


lsvv<-  £o>l  - 

VH 

HV 

HH 


iff) 


I*  ls„i>! 
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.  0' 
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(l-  2531a  +  “316H1  -  <&  -  Oiozloz]  +  t«31c  -  Kz  +  'W^lc*  ) 


where 


+  ^Hxiv^z 
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0 

-1 


aoz^%u~  yozHni 


an a  =  2(^-  D 

lift 

21c 

31a 

31ft 

31c 


1 

1 

0 

0 

('Y  +  Ot  ) 

'•oz  oz} 


(IV.24.a) 


“  2ana  +a116  =-  2(v2-  l),and 
2a3l0  +  ^3lft“2^2-  1>fro*+ao*>- 


(IV.24.b) 

(IV.24.C) 


We  can  arrive  at  the  result  (5  1)  by  substituting  (IV.24.a,b,c)  into  (IV.23)  or 
alternatively  by  using  (IV.  16)'  in  (IV. 17)  to  obtain  the  a  elements  directly  with 

2  co 

ft  - - a  ,B  =  a 

ri  — >■  o  *  — >  o 
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and  substituting  those  results  directly  into  (IV. 22)  in  place  of  the  elements  Bjk.  If  we  choose 
the  latter  course,  then  the  resulting  elements  are: 


2gj  ,  „  2oj  „ 

Bjk  =  V  =  'A7T1  “7 

'  *4  O' 


where 


2o>  .  2oj  „ 

-  “  “o)  =  "  ~  1)Q!oxaoy* 


C 

2oj  A 


2oj 


fl2i(  v  ~«o)  =  ~  iK*. 


~  .  2cj  2cj  0  „ 

a  3i(«0;  -  v  “o)  =  “7 (l>  1)[("oz  +  7oz)a°> "  y°*]' and 


2cj 


2a> 


a4l<2o;- 


Substitution  of  (IV.25)  into  (IV.22)  or  the  alternative  approach  through  (IV .23)  and 
(IV.24.a,b,c)  yields  (5.1). 
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